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ABSTRACT 

We present Hi 21 cm observations of the Orion Nebula, obtained with the Karl G. Jansky Very 
Large Array, at an angular resolution of 7'.'2 x 5'.'7 and a velocity resolution of 0.77 km s _1 . Our data 
reveal H I absorption towards the radio continuum of the H n region, and H I emission arising from the 
Orion Bar photon-dominated region (PDR) and from the Orion-KL outflow. In the Orion Bar PDR, 
the H I signal peaks in the same layer as the H2 near-infrared vibrational line emission, in agreement 
with models of the photodissociation of H2. The gas temperature in this region is approximately 
540 K, and the H 1 abundance in the interclump gas in the PDR is 5-10% of the available hydrogen 
nuclei. Most of the gas in this region therefore remains molecular. Mechanical feedback on the Veil 
manifests itself through the interaction of ionized flow systems in the Orion Nebula, in particular the 
Herbig-Haro object HH 202, with the Veil. These interactions give rise to prominent blueward velocity 
shifts of the gas in the Veil. The unambiguous evidence for interaction of this flow system with the 
Veil shows that the distance between the Veil and the Trapezium stars needs to be revised downwards 
to about 0.4 pc. The depth of the ionized cavity is about 0.7 pc, which is much smaller than the depth 
and the lateral extent of the Veil. Our results reaffirm the blister model for the M42 H 11 region, while 
also revealing its relation to the neutral environment on a larger scale. 

Subject headings: Hn regions; ISM: individual (Orion Nebula, NGC 1976, M42, Orion A, Orion Bar) 



1. INTRODUCTION 

The Orion Nebula (M42, NGC 1976, Orion A) is the 
nearest region of recent massive star formation, contain- 
ing the densest nearby cluster of OB stars. Since the 
optically visible nebula M42 is located in front of the par- 
ent molecular cloud OMC-1, it is accessible for detailed 
studies in every region of the electromagnetic spectrum. 
As a result, the Orion Nebula has become a cornerstone 
for our understanding of massive star formation, as well 
as its feedback effects on the star forming environment, 
which is the subject of the present paper. 

The Orion nebula and OMC-1 are located near the 
center of a prominent north-south ridge of dense molec- 
ular gas, shaped approximately like an integral sign 
(iBallv et al.lll987t iCastets et al.lTT990l: iHever et al.l ll992t 
I Johnstone fe Ballvi 119991 : iPlume et al.1 l2000f ). and con- 
taining the OMC-1 through OMC-4 molecular clumps. 
OMC-1 is the most most prominen t of these, with a 
mass of approximately 2200 M (jBallv et al.l I1987D . 
The integral-shaped ridge is the northern part of 
the larger Orion A giant molecular cloud (GMC), 
which has a mass of about 10 5 M (jMaddalena et al.l 
1986) and is one of a system of two GMCs (the 
Orion A and Orion B GMCs, named after the ra- 
dio sources they contain) that extends roughly north- 
south throug h the belt and sword regions of the Orion 
constellation (iKutner et aLlll977 : Maddalena et al.lll986l : 
iSakamoto et al.lll994t iWilson et al.ll2005l ). These clouds 
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are associated with even large r diffuse Hi clouds 
(jChromev et al.lll989t lGreenlll991| ). An excellent recent 
review of star formation and molecular clouds in th e 
greater Orion region has been presented by I Ballvi (2008). 

The Orion A molecul ar cloud hosts s everal generations 
of OB star formation (|Blaauwl 11964f ) , the youngest of 
which is the Orion Nebula Cluster (ONC ), ionizing the 
M42 Hn region (see iMuench et all 12008. for a detailed 
recent review). This cluster has a central density of 
about 2 x 10 4 stars pc -3 and a total stellar mass of about 
1800 Mq in about 3500 stars (jHillenbrand fe Hartmann 
1998). out to a radius of ~ 2.5 pc. The total mass of 
the ONC is therefore comparable to the molecular gas 
mass o f OMC-1, which is 2200 M©, within a similar 
radius ([Bally et al.l I1987D . Locally, the star formation 
efficiency (here quantified as M stars /(M stars + M gas )) is 
therefore quite high at approximately 50%. On the scale 
of the integral-shaped ridge (linear size about 9 pc) , 
which has a gas mass of - 5000 M (jBallv et al.lll987ji . 
this efficiency is somewhat lower, approximately 25%. 
The ionizing luminosity of the ONC is dominated by 
9 1 C Ori. This star is the most luminous component 
of the asterism formed by the Trapezium stars {9 l A- 
D Ori). 9 1 C Ori is an oblique magnetic rotator with 
an effective temperature T B ff ~ 39 000 ± 1000 K and 
logg = 4.1 (jSimon-Diaz et all I2006T) . implying a spec - 
tral type 06Vp. Observations bv iWeigelt et al.l (|1999l ) 
revealed that 9 1 C Ori is a close binary, dominated in 
mass and luminosity by the star 9 Cj Ori, for which a 
spectral type 05.5 was derived bv iKraus et al.1 (|2007f ). 
The ionizing photon flux corresponding to spectral types 
06 t o 05.5 is Qq = 1.0 - 1.3 • 10 49 s" 1 (IMartins et all 
120051) . 

Most visual studies of the Orion Nebula have concen- 
trated on the ~ 5' diameter optically bright region cen- 
tered on the Trapezium stars, commonly referred to as 
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theHuygens region, after its first description bv lHuvgensI 
(1659). However, lower surface brightness nebular emis- 
sion extends significantly towards the southwest. Includ- 
ing this fainter region the nebula subtends an approx- 
imately circular region on the sk y, with a diameter o f 
about half a degree (e.g., Fig. 1 in lMuench et alJ l2008). 
This region is now referred to as the Extended Orion 
Nebula (EON, iGiidel et al.ll2008[) and contains the Huy- 
gens region at its north-east boundary. The Huygens re- 
gion itself is bou nded at the northeast s ide by the North- 
east Dark Lane (jO'Dell fc Harrisll2010t l. which separates 
M42 from the fainter H II region M43 towards the north- 
east. Another prominent dark feature, already seen by 
iHuveend (|1659ft is the Dark Bay, which is a tongue of 
obscuration, covering part of M42 east of the Trapezium 
stars. 

The Orion Nebula is a blister- type H II region, with the 
ionized gas stream ing away from the high pressure inter - 
face with OMC-1 (IZuckermanlll97l IBalick et alJU97l . 
Velocities with respect to the Local Standard of Rest 
(LSR) in the ionized gas are ulsr — 7.4 ± 1.5 km s _1 
for the low ionization lines ([Oi], [S il] ) arising at the 
ionization front (IF), but lower (i.e., more blueshifted) 
velocities ulsr. = — 0.2 ± 1.3 km s _1 are found for 
higher excitation species ([On], [Oni], [Nil]) and for 
the bulk i onized gas t raced by hydrogen r ecombina- 
tion l ines dKaleil [19671: lO'Dell fc Wen] T1 992; Do i et al 
20041: IHennev et all 120071: iGarci'a-Diaz fc HennevI 1200 
Garcia-Diaz et aLll2008f ) . The background molecular gas 
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is at ulsr ~ 10 km s~ 1 ()Lorenlll979r FI. The ionization 
front (IF) separating M42 and OMC-1 is located be- 
hind the Trapezium st ars, at a distance of approximately 
0.3 pc behind 9 1 C Ori (|Wen fc O'Del]||1995l;lO'l5ell205ll: 
lO'Dell et all 120081 ). A three-dim ensional model of th e 
ionized region has been derived bv lWen fc O'Delll ([19951 ). 
who showed that the IF, which is approximately face-on 
in the region behind the Trapezium stars, curves to an 
orientation that is almost edge-on approximately 100" 
south-east of the Trapezium stars. In this region the IF 
is observed as a prominent, almost linear optical feature 
commonly referred to as the Bright Bar. On the molec- 
ular side of the IF, a photon-dominated region (PDR) 
has formed, which is close to edge-on south-east of the 
Bright Bar. It has been studied in the stron g neutral gas 
coolin g lines, in partic ular [C ill 158 urn (|Stacev et al.l 
11993ft and [Oi] 63 /jm (|Herrmann et al.lll997l ) as well as 
numerous other species. Due to its aspect and proximity, 
the edge-on Orion Bar PDR has become the most iconic 
region of its type. 

The OMC-1 molecular cloud behind M42 harbors 
an obscured region of young massive star formation, 
exhibiting luminous infrared emission with a bolo- 
metric luminosity of about 8 • 10 4 L Q (jGezari et al.l 
19981). known as the Kle inmann-Low region (Orion-KL, 



Kleinmann fc Lowlll967l) . and located about V northwest 
of the Trapezium stars. This region contains a complex 
system of outflows and masers, various young stellar ob- 
jects, and the eponymous Orion Hot Core, a compact 
region of molecular gas and dust with high temperature 
(several 100 K) and density (~ 10 6 cm -3 ) driving a com- 
plex chemistry. The high velocity outflow originating in 

5 In the region under consideration in this paper, LSR and he- 
liocentric velocities are related by i>lsr = ^hol — 18-1 km s — 1 . 



this regio n gives rise to the fam ous "fingers", first dis- 
cussed by lAllen fc Burton! (|1993ft . All of these features 
are the subject of a vast literature, to which we will re- 
turn in Sect. 16.31 E xtensive background can be found in 
the review bv IGenzel fc Stutzk i (1989), which contains 
an over view and syn t hesis of earlier results, and the re- 
view bv 10 'Dell et all (12001 . which discusses more recent 
results on this complex region. 

A second active star forming region is located about 1.'5 
south of Orion-KL. This region, referred to as Orion- 
S, has an infrared luminosity of about 10% of that of 
Orion-KL ([Mezger et al.lll990D . Like Orion-KL, Orion-S 
is a rich source of molecular line em ission, containing sev- 
eral hot cores ([Zapata et al.ll2007h and multiple bipolar 
outflows and maser systems. However, unlike Orion-KL, 
Orion-S is an isolated molecular core located within the 

As a re- 



cavity containing the ONC ([O'Dell et al.l l 
suit several of the outflows originating from Orion-S pro- 
duce optically detectable features, ma ny of which are cat- 
alogued as Herbig-Haro (HH) objects (lO'Dell et alj|1997t 
| Bally et al.l [20001 IHennev et alJl200"l lO'Dell fc HeTmevI 
2008). 

In front of the ionized nebula, several layers of pre- 
dominantly neutral atomic gas are found. These were 
first detected in 21 cm Hi absorption towards the neb- 
ular radio continuum (iMullerl 119591: [Clark et al.l 119621: 
ClarUl965HRadhakrishnan et al.lll972tlLockhart fc Goss 



1978 ) , and are collectively referred to as the Veil ([O'Dell 
2001). The term Veil is appropriate since this fea- 



ture is largely transparent, and only becomes opaque 
(at visual wavelengths) in the Dark Bay and Northeast 
Dark Lane regions, where it s column density is highest 
(jO'Dell fc Yusef-Zadeh|[2000l) . The first full Hi aperture 
synthe sis observations of the O rion Nebula were carried 
out bv lLockhart fc Gossl ([1978D at an angular resolution 
of 2', using the Owens Valley Interferometer. These au- 
thors first showed the presence of 3 velocity components 
in the Veil. This velocity structure was confirmed in 
higher resolution (1 6") aperture synthesis using the VLA 
in C-configuration (|Van der Werf fc G oss 1989, hereafter 
vdWG89), who found LSR velocities of approximately 6, 
4, and —2 km s _1 for the absorbing components A, B and 
C (adopting the notation of vdWG89, which we follow 
in the present paper). These observations confirmed the 
physical associati on of the Veil with the Orion Nebula, 
first suggested bv lLockhart fc Gossl fl978), based on the 
increasing H I column density towards the Dark Bay and 
Northeast Dark Lane in the velocity components A and 
B. Absorption by components A and B is also detected 
towards the smaller Hn region M43 towards the north- 
east, confirming that the Veil represents an extended 
layer covering the M42 /M43 system. The total H I opac- 
ity distribution of components A and B correlates well 
with the optical extinction towards the Huygen s region 
(lO'Dell et al.H1992HO'Dell fc Yusef-Zadehl[2000l) . Physi- 
cal co nditions in the Vei l have been studied further by op- 
tical (lO'Dell et alj|1993ft and ultraviolet (UV) ab sorption 
lines (|Abel et al 112004 l200l iLvkins et al.ll2010D . Model- 
ing of these results has resulted in a location for the Veil 
of a significant, but not accurately determine d, distance 
of 1 — 3 pc in front of the Trapezium stars (jAbel et all 
12004ft . 

Several additional H 1 absorption components have 
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been detected towards the Huygens region. These cover 
only small parts of M42 and are not detected towards 
M43. Velocity com ponent C at ut.sr ~ —2 km s -1 was 



already detected by lLockhart fc Gossl (|1978f l . The H I ob- 
servations described by vdWG89 unexpectedly revealed a 
remarkable set of small-scale (0.02 — 0.06 pc) H I absorp- 
tion components (|Van der Werf fc Gossl 119901 hereafter 
vdWG90). Most of these features (D-G in the nota- 
tion of vdWG90) are blueshifted with respect to both 
the molecular and the ionized gas, and have central LSR 
velocities from —7 to —17 km s -1 . Two of the features 
exhibited several velocity components. In addition, one 
feature (H) was detected by vdWG90 in absorption at the 
velocity of the background molecular cloud OMC-1. The 
features are most likely associated with M42 (vdWG90), 
but their precise nature remained somewhat unclear. 

In the present paper, we return to the Orion Nebula 
to investigate the radiative and mechanical feedback of 
the ONC, the Orion Nebula and the various outflow sys- 
tems, on the neutral environment of the nebula. We use 
new high resolution Hi radio observations to probe Hi 
emission from behind the Huygens region and from the 
Orion Bar PDR, as well as Hi absorption from the Veil 
and the small-scale absorption components. We thus ob- 
tain a comprehensive picture of the radiative and me- 
chanical feedback effects of massive star formation in 
this region on the neutral gas environment. We de- 
scribe the observations and data reduction in Sect. [21 
The radio continuum, H I emission, and H I absorption 
results are presented in Sects. EHSl These results are 
discussed in detail in Sects. M an d El Finally, our con- 
clusions are summarized in Sect. [51 Throughout this 
paper, we adopt a dist ance to the Trapezium stars of 
436 ± 20 pc as given by lO'Dell fc Hennevl (|2008f l. which 
is based on a wei ghted combination of several paral- 
lax measurements dGenzel et al.lll98ll:lKraus et al.ll20"07t 



iMenten et aT]|2007t lHirota et alj|2007l) . At this distance 
1" = 0.0023 pc or 436 AU. Where we use distance- 
dependent quantities from earlier publications, we have 
tacitly converted these to the distance adopted here. 

2. OBSERVATIONS AND REDUCTION 
2.1. Observations 

We used the NRAO Karl G. Jansky Very Large Array 
(VLA), to obtain Hi data of the Orion Nebula in two 
periods in 2006 and 2007 (programs AG738 and AV297). 
The C and then the B array of the VLA were used to ex- 
tend the angular resolution, sensitivity and velocity cov- 
erage of the old C array data of vdWG89 and vdWG90, 
obtained in 1984. The VLA correlator was used. The to- 
tal bandwidth was 781.25 kHz with 256 channels and two 
circular polarizations, centered at ^lsr = 2.0 km s . 
The channel separation was 3.052 kHz (0.64 km s -1 at 
the Hi line) and the velocity resolution was 0.77 km s . 

The C array data was obtained in a series of three 
5 hour observations on 2006 September 29, November 9 
and November 19. Twenty of the VLA antennas were 
used, with no use of the 7 antennas that had been con- 
verted to the EVLA at that time. The phase calibration 
was based on frequent observations (once per half hour) 
of the quasar OG 050 ( J0532+0732) with a flux density 
of 1.8 Jy. The flux density scale was set by observations 
of 3C48 (15.7 Jy). 



The B array data were obtained during late 2007 in a 
three times 5 hour observation on November 15, Novem- 
ber 24 and December 3. The flux density scale was set 
using observations of 3C 138 (total flux density 8.3 Jy). 
The observations of 3C 138 were carried out every 30 min- 
utes for a period of 4 minutes. 

For both the C array and the B array data, bandpass 
responses of each antenna were determined by observing 
the strong sources 3C48 and 3C84; these observations 
were shifted by plus and minus 0.7 MHz (148 km s -1 ) to 
avoid the H I emission near ^lsr = km s _1 and absorp- 
tion lines of Galactic H I in the spectra of the calibration 
sources. 

2.2. Reduction and generation of data cubes 

During the 2007 observations, we used EVLA anten- 
nas for the first time. At this time there were 12 EVLA 
antennas and 13 VLA antennas. During a test obser- 
vation of 3C48 obtained on 2007 October 4, an alias- 
ing problem with the old VLA correlator and the use of 
EVLA baselines was discovered by a number of NRAO 
staff (including M. Goss). This problem was caused by 
the hardware used to convert the digital signals from the 
EVLA antennas into analogue signals to be fed in the 
VLA correlator, which caused power to be aliased into 
the bottom 0.5 MHz of the baseband. Only EVLA to 
EVLA antenna correlations were affected. A number of 
partial solutions were founc0. For the Orion A H I obser- 
vations, the solution adopted was to use observations of 
the strong source 3C 138 every 30 minutes and to apply a 
time variable baseline based calibration (as opposed to an 
antenna based calibration) to correct for the closure er- 
rors due to the mismatched and time variable bandpasses 
resulting from the aliasing. This scheme was tested in de- 
tail using the correlator configuration that we used for 
the Orion A Hi observations. We found that tracking 
the errors over this time interval worked well and the 
visibility functions of the calibrator sources for EVLA to 
EVLA baselines had a similar behavior as those of the 
VLA to VLA baselines. Before the phase closure cor- 
rections were made, the amplitude fluctuations were at 
the level of 15% (la) of the continuum flux density of 
3C 138; after correction the fluctuations were reduced to 
values well below 2%. With the advent of the WIDAR 
correlator in early 2010, these aliasing problems have dis- 
appeared. 

The data from the C and B arrays were then combined 
and the line images were made after subtracting the con- 
tinuum in the uv plane (using the AIPS task UVLSF); 
102 of the 255 channels were line free and formed the 
continuum. The final images were made using the AIPS 
task IMAGR with Robust=0 weighting. The resulting 
datacube has a synthesized beam of 7'/2 x 5'/7 at a po- 
sition angle of 29.7° and an r.m.s. noise per channel of 
40 K. The conversion factor between brightness tem- 
perature and flux density is 14.9 K (mJy beam -1 ) -1 . 
A 1420.4 MHz image was produced using a multi- 
scale CLEAN algorithm, in order to optimally preserve 
the large range of spatial scales present in this image. 
The measured r.m.s. noise in the continuum image is 
3.1 mJy beam -1 . 

6 http:/ /www. via. nrao.edu/astro/guidcs/cvlareturn/aliasing/ 
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In the spectral line data cube, we discarded channels 
with elevated noise at the edges of the band. Our final 
data cube covers the range —62.00 km s _1 < ulsr < 
68.56 km s _1 . 

2.3. Further processing 

Inspection of the Hi data cube revealed a large and 
complex set of features at various velocities, and with 
various angular sizes. Remarkably, H I is detected in both 
emission and absorption. 

Interferometric imaging of extended low-level emission 
features in the presence of a strong continuum requires 
careful processing, because of non-linearities introduced 
by deconvolution algorithm s such as CLEAN (see e.g., 
I Van Gorkom fc EkerslU989l ) . which may give rise to spu- 
rious features after continuum subtraction. The best 
way to avoid these problems consists of first subtract- 
ing the continuum and applying the deconvolution to the 
continuum-free line images. As described above, this is 
the procedure that was used. After continuum subtrac- 
tion the H I absorption produces a strong negative signal 
carrying the imprint of the subtracted continuum at LSR 
velocities between —20 and 10 km s _1 ; at other veloci- 
ties any remaining signal results from Hi emission. As 
a result, Hi emission can only reliably be studied at ve- 
locities outside the range from —20 to 10 km s _1 . In 
order to increase the S /N ratio of the H I emission data, 
we convolved the channel maps to a 7'.'5 circular beam, 
and smoothed the data cube spectrally by a factor 2, 
i.e., to a velocity resolution and channel separation of 
1.29 km s . The r.m.s. noise in brightness temperature 
lb in these images is 20 K. Since the shortest baselines 
in our observations were about 73 m, our data are insen- 
sitive to structures with scales of about 7' or more. 

In order to study the H I absorption, the full spatial and 
spectral resolution H I line data cube was used, with the 
corresponding 1420.4 MHz continuum image, to derive a 
data cube of H I optical depth r, following the approach 
of vdWG89. Optical depths were derived by solving the 
equation of transfer 



T b («) = [1 - e- T ^][T s -T c - T back (t>)], 



(1) 



where T\>(v) is the observed Hi brightness temperature 
at LSR velocity v, after subtraction of the continuum 
(which has brightness temperature T c ). T s is the spin 
temperature of the absorbing Hi, and Tb ac k(v) is the 
brightness temperature (at LSR velocity v) of Galac- 
tic Hi originating behind the absorbing Hi. The peak 
brightness temperature of Galactic Hi in the re gion of 
the Orion Nebula is about 60 K (jGreenl H99ll ). It is 
not possible to determine what fraction of this signal 
originates behind the absorbing H I, and the situation is 
complicated further by the fact that this fraction may 
be a function of v. Therefore the observed Galactic Hi 
brightness temperature only provides an upper limit for 
?back(v). For T s a harmonic mean value can be deter- 
mined at positions where Hi 21 cm absorption can be 
combined with measurements of Lya absorption. Such 
measurements are avail able at the pos itions of 9 1 C Or i 
UShuping fc SnowlH997l) and 1 B Ori (I Abel et al.H200l) . 
giving T s » 90 K (80-110 K) in component A and 
T s w 135 K (100-160 K) in component B. Given the 
uncertainties in T s and Tb ac k(w) we solve Eq. ([T]) with the 



approximation T c 3> \T S — Tback(^)|, by only calculating 
opacities at positions where T c > 450 K (correspond- 
ing to a surface brightness level of 30 mJy beam - ), 
which is almost 10c in the continuum image. While Hi 
absorption lines can be detected towards considerably 
fainter continuum levels, quantitatively reliable opacities 
can only be derived where T c J> \T S — Tb ac k(u)|. With 
this approach the precise value of T s for determining the 
opacities (as well as the implicit assumption of Eq. (fTJ) 
that T s is the same at every position and for all absorbing 
velocity components) becomes irrelevant. At positions 
with T c < 450 K, this procedure will give rise to sys- 
tematic errors in the derived opacities; at these positions 
opacities are therefore not calculated. Q 

3. CONTINUUM EMISSION 

The 1420.4 MHz continuum image is shown in Fig. [TJ 
This image shows both the main Hn region M42 and 
the fainter Hn region M43 in the northeast. The im- 
age of M42 is dominated by the Huygens region and 
the Bright Bar in the southeast. Fainter emission from 
the EON is seen to extend in all directions from the 
southeast counterclockwise to the West, but not in the 
other directions. Indications for the presence of low 
surface bright ness radio emission f rom the EON were 
first fo und bv lMills fc Shaven (11968ft and lGoss fc Shaverl 
(1970). Emission from the ful l EON is detected i n the 
4.75 GHz single-dish image of IWilson et al.l (11997f) and 
in the 330 MHz VLA image of ISubrahmanvan et al.l 
(12001; the 1.5 GHz VLA imag e by lYusef-Zadehl (|1990t ) 
and ISubrahmanvan et al.l (|2001l ) also shows the extended 
emission from the EON. 

The total flux density of M42 in our image is 335 ± 
15 Jy. This value is somewhat lower than the total flux 
density of 374 Jy found by vdWG89 (which is consistent 
with the best single-dish value, see Table 2 in vdWG89). 
For M43, we find a total flux density of 14 ± 2 Jy, in 
excellent agreement with vdWG89. 

The peak continuum flux density of M42 is 
376 mJy beam -1 , which corresponds to a peak contin- 
uum brightness temperature T c = 5600 K. The electron 
temperature in this region is T e = 8400 ± 400 K, as de- 
termined by IWilson et al.l (|1997l ) from measurements of 
the H64a recombination line. Since T c = T c (l — e _Tff ), 
the peak free-free optical depth is rg = 1.1, i.e., the 
Hn region is significantly optically thick at 1420.4 MHz, 
and will be opaque at lower frequencies, in agreement 
with the spectral index d istribution between 330 an d 
1500 MHz determined by ISubrahmanvan et~aH (|200lD . 
A free-free opacity rg w 1 is also found at the brightest 
peaks of the Bright Bar. 

4. Hi EMISSION FEATURES 

4.1. Hi emission images 

As noted in Sect. 12.31 Hi emission can be studied at 
velocities avoiding strong absorption, i.e., outside the ve- 
locity interval from —20 to 10 km s . Since the promi- 
nent molecular cloud associated with the Orion Nebula 

7 The derived H I opacity datacube and the corresponding 21 cm 
continuum image will be made available in electronic form through 
the Centre de D onnees astronomiques de Strasbourg (CDS) at 
http : //cds . u-strasbg . f r 
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50 3 40 s 30 3 30 s 5 h 35 m 10 s s 34 m 50 s 

E.A. [J2000.0] 

Fig. 1. — Continuum emission of the Orion Nebula at 1420.4 MHz, constructed using a multi-scale deconvolution (see Sect. 12. 2"t . Contours 
indicate surface brightness levels of 20, 40, 60, 80 and 100 mjy beam - 1 (black contours) and 150, 200, 250, 300 and 350 mjy beam -1 
(white contours). The small ellipse in the lower lefthand corner indicates the FWHM size and the orientation of the synthesized beam 
(7'.'2 X 5'.'7 at a position angle of 29. 7% The image has been corrected for primary beam attenuation. The principal massive young stars 
are indicated by red dots. A cyan dot indicates the position of the Orion-KL region. 



is at fLSR ~ 10 km s -1 , Hi emission from the neutral 
environment of the Hil region can be probed only in 
its red line wing. Inspection of the Hi emission data 
cube revealed Hi emission at LSR velocities from 10 to 
19 km s _1 , and six images covering this velocity range 
are presented in Fig. [3] 

Inspection of Fig. [2] reveals that the strongest Hi emis- 
sion is found in the region directly southeast of the Bright 
Bar. H I brightness temperatures in this region are ap- 
proximately 250 K (coded red in Fig. [5]), with peaks 
reaching 300 K, indicating that the gas is quite warm. 
While the brightest Hi emission is found closest to the 
Bar, the emission extends towards the southeast over a 
distance of about 6' (0.8 pc), at brightness temperature 
levels of about 120 K. Other features worthy of note are 
an elongated H I feature extending from the Bar region 
towards the southwest at velocities of 13 — 16 km s^ 1 , 



and compact H I emission features with higher velocities 
(^lsr up to 19 km s _1 ). Hi emission features northeast 
of M42 trace the direct environment of M43. 

4.2. Position-velocity diagrams of Hi emission 

In order to study the velocity structure of the H I emis- 
sion and its relation to the Hi absorption (to be dis- 
cussed below), we have constructed a number of position- 
velocity (PV) diagrams of the H I emission. The orienta- 
tions of the spatial axes of these diagrams are shown in 
Fig. El 

Figure [4] shows the velocity structure of the extended 
H I layers in the region of the Orion Nebula. A prominent 
H I layer can be seen in emission in the south (at an offset 
of approximately —300" in Fig. [3}. Following this layer 
northwards, it produces strong Hi absorption in front of 
the strong radio continuum of M42. Between M42 and 
M43, in the Northeast Dark Lane, strong Hi emission 
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50 s 30 s 5 h 35 m 10^34 m 50 3 50 s 30 s 5 h 35 m 10^34 m 50 s 

R.A. [J2000] 

Fig. 2. — Images of Hi emission at the LSR velocities indicated in the upper left-hand corner of each frame. The images cover a velocity 
interval of 1.28 km s — 1 each. Colors ranging from dark blue to red indicate a brightness temperature range from to 250 K. White contours 
indicate the 21 cm continuum at levels of 50 and 200 mjy beam -1 . These images have been corrected for primary beam attenuation. The 
7'.'5 circular beam of these images is indicated in the lower righthand corner. The red dots in the lower left panel indicate the principal 
young massive stars (see Fig. [T] for legend). 
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Fig. 3. — Location of the position- velocity (PV) diagrams of Hi emission shown in Figs, l4lf71 superposed on an Hi emission image at 
^LSR = 14.9 km s — 1 (greyscale). The precise lengths, orientations and positions of these PV diagrams are indicated by the colored bars. 
In each bar, a thick dot of the same color indicates the zero position of the spatial coordinate in the corresponding PV diagram. The 
relevant figure numbers are indicated. White contours indicate 1420.4 MHz continuum surface brightnesses of 15, 40, 70, 100, 150, 200 and 
300 mjy beam -1 . The purple rectangle in the Bright Bar region indicates the region used for constructing the s trip shown in Fig. 1191 The 
yellow rectangle (which is contained in the purple one) shows the region over which the spectrum shown in Fig. 1181 was averaged. 
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offset ["] 

Fig. 4. — Position-velocity diagram of Hi emission along the red 
north-south bar shown in Fig. \3\ The spatial axis has its zero posi- 
tion at RA = 5 h 35 m 29?5, Dec = -5° 23' 37" and a position angle 
of 0° Spatial offsets are positive towards the north and negative 
towards the south. Features discussed in the text are indicated. 




NE -300 -200 -100 100 200 300 STf 

offset ["] 

Fig. 5. — Position- velocity diagram of H I emission along the light 
blue northeast-southwest bar shown in Fig. [3] The spatial axis has 
its zero position at RA = 5 h 35 m 03?5, Dec = -5° 27' 40" and a 
position angle of 64.66° Spatial offsets are negative towards the 
northcst and positive towards the southwest. Features discussed 
in the text are indicated. 
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offset ["] 

Fig. 6. — Position- velocity diagram of H I emission along the dark 
blue northwest-southeast bar shown in Fig. [3] The spatial axis has 
its zero position at RA = 5 h 35 m 23?8, Dec = -5° 24' 38" and a 
position angle of —45° Spatial offsets are negative towards the 
southeast and positive towards the northwest. Features discussed 
in the text are indicated. 



is found. These emission features have central velocities 
VhSR ~ 2 km s -1 . The absorbing Hi in front of M43 is 
at ^lsr ~ 7 km s _1 . 

In the region southwest of the Bright Bar (offsets —160 
to 0" in Fig. [4]) H I emission is found with a peak velocity 
of about 12 — 13 km s _1 , i.e., displaced in velocity from 
the Hi absorption by about 10 — 11 km s _1 . This Hi 
emission feature is also detected at offsets —300 to —100" 
in Fig.[5j which presents a PV diagram through the elon- 
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offset ["] 

Fig. 7. — Position- velocity diagram of Hi emission along the 
green northwest-southeast bar shown in Fig. [3] The spatial axis 
has its zero position at RA = 5 h 35 m 25?4, Dec = -5° 22' 48" and 
a position angle of —80° Spatial offsets are negative towards the 
southeast and positive towards the northwest. Features discussed 
in the text are indicated. 

gated Hi feature detected at 14 km s _1 in Fig. [2] This 
PV diagram clearly reveals the elongated Hi feature as 
a kinematically separate entity, detected at offsets from 
280 to —80". At the latter position, it connects to the 
Hi at 12 km s _1 southeast of the Bright Bar. 

Figure [5] shows a PV diagram crossing the Orion Bar 
PDR orthogonally, with the IF at offset 0", and also 
crossing the compact high- velocity H I emission feature 
detected in Fig. [2] at «lsr ~ 18 km s _1 . The Hi emis- 
sion at «lsr ~ 14 km s _1 is located in the region of the 
Orion Bar PDR (offsets -220 to 0"), but also extends 
slightly northwest of the IF (offsets to 50"). The fea- 
ture M is an absorption feature associated with the IF 
that will be discussed in Sects. 15.11 and 15.3.91 The high 
velocity H I feature at offsets 100 — 120" reaches velocities 
up to 31 km s^ 1 , and contains two distinct components 
separated by a local emission minimum. This structure 
is well detected in Fig. [JJ In this diagram, faint high 
velocity H I emission from the northwest component (off- 
set 100") is also detected at negative velocities. The PV 
diagram in Fig. [Jj also crosses a region of H I emission at 
15 km s" 1 (at offsets from —30 to 30") located at the 
position of the optical Dark Bay. This feature is remark- 
able since it contains at its center an absorption feature 
(marked I in Fig. [Jj and discussed further in Sect. 15.3.9)) . 
Figure [7J also shows Hi emission at Ulsr ~ 12 km s _1 at 
offsets from -160" to 120", associated with OMC-1, but 
strong foreground H I absorption precludes further study 
of this feature. 

5. H i ABSORPTION FEATURES 

5.1. Overall velocity structure of the absorbing Hi 

The velocity structure of the Hi absorption is illus- 
trated in Fig. [8l which shows the H I opacity spectrum av- 
eraged over a region in the southeast part of the Huygens 
region, close to the Bright Bar. This spectrum matches 
that shown in Fig. la of vdWG90, corresponding to ap- 
proximately the same region. Only unsaturated points 
were included in the calculation of the average spectrum. 
Therefore, the spectral shape in the vicinity of the peaks 
of components A and B in this figure should be treated 
with caution. 

Five prominent H I velocity components are seen in this 
figure, corresponding to the components A, B, C, Dl and 
D2 of vdWG89 and vdWG90. The H I components A and 
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Fig. 8. — Spectrum of Hi opacity, averaged over a 27" X 27" area 
centered on RA=5 h 35 m 22?88, Dec=-5° 24' 30'.'7, as indicated in 
Fig. [9] in the panel at dlsr = —22.5 km s _1 . 

B cover the entire nebula. The difference in the central 
velocities of these two components decreases towards the 
northeast. As a result, and due to the increasing opacity 
towards this region (resulting in saturation of the line 
peaks), components A and B become difficult to sep- 
arate in the region towards the Northeast Dark Lane. 
Over most of the nebula, these components can however 
be traced as two kinematically distinct features. Compo- 
nent C is prominent in the southwest part of the Huygens 
region. It is however not detected in the region of the 
Trapezium stars and towards the Northeast Dark Lane. 
Towards M43 (even farther to the northeast) only the 
components A and B are detected. 

Figure E] also shows two velocity components, Dl and 
D2, at significantly negative velocities. These are ex- 
amples of features with small spatial scales compared to 
components A and B, as discussed by vdWG90. Two of 
these features (D and F) revealed several velocity com- 
ponents. Because of the limited frequency coverage of 
the observations by vdWG90, components D2 and F2 
were not observed over their full velocity extent. Fig- 
ure [5] shows that D2 is completely within the spectral 
band of the observations presented here; this is also the 
case for F2 (not shown in Fig. [8]). In addition, our obser- 
vations enable the detection of further Hi components 
at velocities not covered by vdWG90. One new veloc- 
ity component is indicated in Fig. [5] as component M. 
This feature, which is also indicated in Fig. [5]) will be 
discussed in Sect. 17.5.21 The detection of this very faint 
feature illustrates the high sensitivity and dynamic range 
of this dataset. 

5.2. Hi opacity images 

In order to study the morphology of the Hi opacity 
as a function of ulsR: we present a sequence of opac- 
ity images. For presentation purposes, the data cube 
was spectrally smoothed (ignoring saturated pixels) to 
a channel resolution and separation of 1.28 km s _1 . A 
sequence of these images is shown in Fig. [9] We note 
that the smoothing in velocity was only done for the 



purpose of creating these figures; all analysis was done 
on the full velocity resolution data cube. The range 
0.6 km s _1 < ulsr < 5.9 km s _1 is heavily affected 
by saturation. Therefore, this region, covering part of 
component A and most of component B, was omitted 
from Fig. HI 

5.2.1. Large-scale Hi absorption features 

We focus first on the components A, B and C identified 
by vdWG89. Referring to Fig. [8j the red line wing of 
component A is traced in the opacity image at i>lsr = 
8.4 km s _1 . A strong increase in opacity is observed 
towards the Northeast Dark Lane, in agreement with the 
results of vdWG89. 

Component C of vdWG89, located at the southwestern 
edge of the Huygens region, can be traced in the opacity 
images between —4.4 and —0.6 km s _1 . The morphol- 
ogy of this component is obvious at «lsr = —3.2 km s _1 , 
where it displays a striking arc of high H I opacity trac- 
ing the extreme southwest edge of the Huygens region, 
and broken into a number of H I peaks. Several other H I 
features are detected in the same velocity range, most no- 
tably near the Bright Bar and towards the northern part 
of Huygens region. However, as will be discussed below 
(Sects. I5".3.5l and l5.3.6[) . those features are not physically 
related to the prominent H I opacity arc at the southwest 
edge of the Huygens region. Therefore in our nomencla- 
ture component C will only denote this Hi arc. 

5.2.2. Small-scale Hi absorption features 

To trace the small-scale Hi absorption components, 
we begin at the most negative velocities and follow the 
opacity images towards more positive velocities. Central 
positions and velocity ranges over which these features 
are detected are summarized in Table [TJ The components 
are also indicated in Figs.l9landlT0l Table [T] also summa- 
rizes Hi masses, and peak and average column densities, 
as well as approximate total sizes of the various features. 
The Hi masses and column densities scale with T s , here 
assumed to be 100 K (see Sect. 17. l|) . For convenience, the 
formulae relating the Hi emission and absorption mea- 
surements for Hi column densities, and the underlying 
assumptions, are summarized in the Appendix. 

The most negative velocity components in our data 
are components D (in the region of the Bright Bar) 
and F (in the western part of the Huygens region). 
Both are detected over a considerable velocity range 
(—21.2 < t>LSR < —10.9 km s _1 ). The more opaque com- 
ponent D displays obvious extended structure with em- 
bedded higher opacity clumps. Both components contain 
several velocity components, in agreement with vdWG90. 

Continuing to more positive velocities, component E 
of vdWG90 is detected in the region of the Dark Bay. 
The data shown by vdWG90 already indicated that this 
component consists of several clumps, which are more 
prominent in the present higher resolution data. Inspec- 
tion of Fig. [9] in the region of component E reveals the 
presence of 10 compact Hi opacity clumps. At the res- 
olution of ~ 6", the diameters of the unresolved clumps 
are less than 0.014 pc or 2800 AU. 

The remaining blueshifted small-scale feature identi- 
fied by vdWG90, component G, is detected as a conspic- 
uous elongated feature towards the northern part of M42 
at velocities from —4 to —11 km s _1 . 
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Fig. 9. — Images of Hi opacity (color images), at the LSR velocities indicated in the upper left-hand corner of each frame. White contours 
indicate the background 21 cm continuum at levels of 50, 100, 200, 300 and 400 mjy beam . Black pixels within the faintest continuum 
contour denote positions affected by saturation. The opacity scale is indicated by the color bar in the first frame. In this frame also the 
region over which the opacity spectrum shown in Fig. [8] has been averaged, is shown. The size and orientation of the synthesized beam 
(7'.'2 X 5'.'7 at a position angle of 29. 7°j is indicated by the ellipse in the top righthand corner of the first frame. Letters C— M indicate 
various absorption features discussed in the text. 
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Fig. 9. — (continued) 
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Fig. 9. — (continued) 
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TABLE 1 

Global parameters of small-scale Hi components 



Name R.A. a Dec. a «LSR range Diameter N(H i) ( j ^Hl ( iocTk ) Associated with 

[J2000] [kms- 1 ] [pc] [10 19 cm" 2 ] [M ] 

peak b mean 



D c,d 


5 h 35 m 23!0 


-5° 24' 


30" 


-21.2 - 


•> -10.9 


lf5 


0.21 


100 


6.6 
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6 2 B Ori 


E c 


5 h 35 m 25?0 


-5° 23' 
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0.02 


HH 202 


Gl f 


5 h 35 m 14?0 


-5° 21' 


40" 


-9.6 - 


■> -1.9 


l.'S 


0.21 


> 350 


3.6 


0.04 




G2 f 


5 h 35 m 13!0 


-5° 21' 


20" 


-12.2 - 


-> -8.3 


30" 


0.07 


> 350 


12 


0.02 


H2 finger 


G3 f 


5 h 35 m 17?5 


-5° 21' 


40" 


-3.2 - 


-> -0.6 


l'.5 


0.21 


27 


2.7 


0.03 




H 


5 h 35 m 13f0 


-5° 24' 


00" 


+7.2 - 


► +12.3 


40" 


0.09 


34 


9.9 


0.02 


Orion-S 


I 


5 h 35 m 25!6 


-5° 22' 


49" 


+9.7 - 


► +13.0 


8" 


0.02 


23 


20 


0.002 


extinction Knot 2 


J f 


5 h 35 m 22?0 


-5° 24' 


52" 


-4.4 - 


* -1.2 


30" 


0.07 


> 340 


3.3 


0.04 


HH 203/204 


K 8 


5 h 35 m 23?0 


-5° 25' 


45" 


-5.7 - 


■> -0.6 


10" 


0.02 








shock south of HH 204 


L d 


5 h 35 m 20?4 


-5° 24' 


25" 
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a Positions for extended features are approximate center positions. 
b Lower limits are affected by saturation. 
c Two velocity components 
d Arc 

e Ten compact clumps 
f Elongated 

s Column density and mass cannot be calculated due to saturation. 
h One of several compact clumps at the edge of the Bright Bar 

One compact feature was identified by vdWG90 at 
positive ulsr (their component H). This feature is at 
the velocity of OMC-1 and it is detected at wlsr = 
7 — 12 km s _1 in Fig.|9]in the western part of the Huygens 
region. 

The data cube was carefully searched for additional ve- 
locity components. No additional H 1 components were 
found at velocities more negative than that of compo- 
nent D, even though Nai and Can absorption line mea- 
surements towards the Trapezium stars and 6 2 A Ori 
reveal several of components at these velocities, e.g., 
at vt.sr ~ ~ 32 - 8 ; -20.6, -17.0, and -12.1 km s" 1 
(|Q'Dell et al.l ll993). Likewise, no absorption was found 
at velocities more positive than that of component H. 
However, our search revealed several new small-scale fea- 
tures within the velocity range shown in Fig. [9] spatially 
and kinematically distinct from the features described 
above. Some of these features are quite compact, or are 
located close to other features, which accounts for their 
non-detection in the lower resolution data of vdWG90. 
Several features are quite close in velocity to compo- 
nent C, but spectra and position-velocity diagrams re- 
veal that these features are kinematically distinct. Their 
global properties are summarized in Table [TJ These fea- 
tures are best seen in PV diagrams discussed below. 

5.3. Position- velocity diagrams of Hi opacity 

The location of a set of illustrative PV diagrams is 
indicated in Fig. 1101 These diagrams are presented in 
Figs. [TTirTTl The various Hi components are indicated 
in these figures. 

5.3.1. Components A, B and C 

Figure [11] shows the saturated signal of the large-scale 
components A and B at positive velocities. In addition, 
component C is detected at slightly negative velocities 



(—4 to —2 km s v ). As already indicated by the channel 
map at —3.2 km s _1 , Fig. ITT1 shows that this component 
does not cover the entire nebula (in contrast to compo- 
nents A and B), but is found only at the edge of the 
Huygens region, where it forms a large arc in the shape 
of an incomplete semicircle. 

5.3.2. Structure and extent of component D 

The velocity structure and extent of the most extended 
small-scale component D is shown in Figs. [T^] and [TBI 
which show crosscuts through this feature in orthogonal 
directions. Figure[12] which shows a PV diagram approx- 
imately along the Bright Bar, shows that this component 
is extended over about 160". The prominent high opac- 
ity features seen at the most negative velocities appear 
to be connected to the main H I components through 
gas at intermediate velocities. Inspection of Fig. [131 
which presents a PV diagram roughly perpendicular to 
the Bright Bar, shows a velocity gradient, in the sense 
that the most negative velocities occur to the southeast. 
The overall velocity structure therefore resembles that of 
an expanding shell. This impression is reinforced by the 
morphology of the absorbing Hi, in particular at veloc- 
ities between —17.3 and —10.9 km s _1 , and also illus- 
trated by the dark blue contours in Fig. [TO] Figure [13] 
also shows that the blueshifted gas of component D con- 
nects to the larger scale Hi layers through gas at inter- 
mediate velocities, located approximately 20" east of the 
Trapezium stars. 

5.3.3. Component E: Hi in the Dark Bay 

A PV diagram through component E, which is located 
in the area of the Dark Bay, is shown in Fig. [TJ] While 
this component is resolved into a number of compact 
opacity peaks, these are obviously part of one coherent 
velocity structure. This diagram also demonstrates that 
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Fig. 10. — Location of the position-velocity (PV) diagrams of Hi optical depth shown in Figs. fTl~Ul7l superposed on an HST/WFPC2 
image (greyscale). The lengths, orientations and positions of these PV diagrams are indicated by the colored bars. In each bar, a thick dot 
of the same color indicates the zero position of th e spatia l coordinate in the corresponding PV diagram. Finally, contour sets indicate the H I 
opacity of the various Hi features shown in Figs. rTTIll7l These contours are taken from the opacity images at —2.5 km s" 1 (yellow contours 
in southwest, indicating component C, and red contours near the Bright Bar, indicating components J and L); —17.3 km s _1 (dark blue 
contours near the Bright Bar, indicating component D, and green contours indicating component F); —7.0 km s — 1 (yellow, component E); 
—7.7 km s _1 (purple contours in northwest, indicating component Gl); —9.6 km s^ 1 (light blue, component G2); —1.3 km s" 1 (red contours 
in the north, indicating component G3); 8.4 km s — 1 (light blue contours southwest of the Trapezium stars, indicating component H); 
10.3 km s — 1 (green, component I); and 1.9 km s — 1 (light blue contours in the south, indicating component K). The optical HST image is 
a combination of the WFPC2 F502N and F658N filter images from lO'Dell fc Wonel 1T996T ), which contain respectively the [Om] 5007 A 
and the [Nil] 6583 A and Ha lines, thus the extremes of both high and low excitation levels. The features indicate d by "SW Cloud" and 
"Knot 1—3" arc features of enhanced extinction identified in the extinction study of O 'Dcll & Yuscf-Zadch (2000). "Dyn-ctr" indicates 
the dyanmical center of the Orion-KL molecular outflow, discussed further in Sect. 16.31 "H2CO absorption cloud" denotes the outline of 
the H2CO absorp tion s ignal of Orion-S. The optical outflow source associated with Orion-S is indicated by "OOS". These features are 
discussed in Sect. 1775,11 
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Fig. 11. — Position- velocity diagram of Hi opti cal depth along 
the yellow southeast-northwest bar shown in Fig. 1101 towards the 
southwest of the Huygens region. The spatial axis has its zero 
position at RA = 5 h 35 m 10?5, Dec = -5° 24' 56" and a position 
angle of 120° Spatial offsets are negative towards the southeast and 
positive towards the northwest. Velocity components discussed in 
the text are indicated. 
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Fig. 14. — Position- velocity diagram of Hi optical depth along 
the yellow sout hea st- northwest bar west of the Trapezium stars, 
shown in Fig. 1101 The spatial axis has its zero position at 
RA = 5 h 35 m 22?4, Dec = -5° 21' 58" and a position angle of 
345°. Spatial offsets are negative towards the southeast and posi- 
tive towards the northwest. Velocity components discussed in the 
text are indicated. 




offset ["] 

Fig. 12. — Position- velocity diagram of H I optica l de pth along the 
dark blue northeast-southwest bar shown in Fig. 1101 The spatial 
axis has its zero position at RA = 5 h 35 m 22?6, Dec = -5° 24' 25" 
and a position angle of 45° Spatial offsets arc negative towards 
the northeast and positive towards the southwest. Velocity com- 
ponents discussed in the text are indicated. 
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Fig. 15. — Position- velocity diagram of Hi opt ical depth along 
the green southeast-northwest bar shown in Fig. 1101 The spatial 
axis has its zero position at RA = 5 h 35 m 13?2, Dec = -5° 23' 04" 
and a position angle of 100° Spatial offsets are negative towards 
the southeast and positive towards the northwest. Velocity com- 
ponents discussed in the text are indicated. 




offset ["] 

Fig. 13. — Position- velocity diagram of Hi optica l de pth along 
the dark blue southeast-northwest bar show n in Fig. 1101 The spa- 
tial axis has the same zero position as for Fig. 1121 but has a posi tion 
angle of 135° perpendicular to the PV diagram shown in Fig. 1121 
Spatial offsets are negative towards the southeast and positive to- 
wards the northwest. Velocity components discussed in the text 
are indicated. 

although in projection components D and E are almost 
contiguous, component E is actually a kinematically sep- 
arate feature. 

5.3.4. Stucture of the high velocity component F 

The PV diagrams of the high velocity component F in 
Figs. [15] and [16] reveal clearly that the high (negative) 
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Fig. 16. — Position- velocity diagram of H I o ptic al depth along the 
green northeast-southwest bar shown in Fig. 1101 The spatial axis 
has the same zero position as for Fig. 1151 but has a po sitio n angle 
of 10° perpendicular to the PV diagram shown in Fig. 1151 Spatial 
offsets are negative towards the northeast and positive towards the 
southwest. Velocity components discussed in the text are indicated. 



velocity gas is connected with the main H I components 
through features at intermediate velocities. This is par- 
ticularly clear in Fig. [TS] where the high velocity gas 
(which shows several velocity components) is connected 
to the lower velocity gas through a prominent H I feature 
at its northwestern side. Careful inspection of Fig. 1151 
and of the Hi opacity data cube reveals that a fainter 
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Fig. 17. — Position- velocity diagram of H I optica l de pth along the 
light blue southeast-northwest bar shown in Fig. 1101 The spatial 
axis has its zero position at RA = 5 h 35 m 14?9, Dec = -5° 22' 01" 
and a position angle of 328° Spatial offsets are negative towards 
the southeast and positive towards the northwest. Velocity com- 
ponents discussed in the text are indicated. 

connection between the high velocity gas and the large- 
scale components is also present at the southeast side of 
the feature. 

The fact that the location of component F spatially 
coincides with a gap in the extended gas layer at less 
negative velocities is striking. This behavior argues that 
component F is physically part of the larger scale Hi 
layer, but that it has been accelerated to negative veloci- 
ties. The situation is further illustrated by Fig. [TBI which 
shows a PV diagram along a position angle perpendicu- 
lar to Fig. [151 The connection with the lower velocity gas 
is clearly observed on both sides of the feature, as well 
as the gap in less negative velocity gas at the position 
component F. 

Like component D, the velocity structure of compo- 
nent F resembles that of an expanding bubble or shell, 
although this interpretation ignores the fact that at some 
positions several velocity components are present. How- 
ever, component F is much more confined than compo- 
nent D, with a diameter of about 60" (0.14 pc). 

5.3.5. The elongated component G 

A PV diagram over the long axis of the elongated 
component G is shown in Fig. 1171 Inspection of the 
line images in Fig. |9] shows that multiple elongated 
H I absorption features are present towards the north- 
ern part of M42 at the full velocity range from —0.6 
to —12.2 km s -1 , which we collectively denote compo- 
nent G. The PV diagram in Fig. [17] shows a prominent 
feature at velocities of —3 to —4 km s _1 . However, at 
more negative velocities (approximately —10 km s _1 ), an 
additional feature is detected towards the northern part 
of the nebula. We denote the features at ~ — 3 km s _1 
Gl and the feature at ~ —10 km s _1 G2, as indicated 
in Fig. [T7] The latter feature matches component G of 
vdWG90. Finally, we note that a parallel and similarly 
elongated Hi absorption feature appears ~ 45" east of 
component Gl at velocities of —1.9 and —0.6 km s™ 1 in 
Fig. [HI we denote the eastern feature component G3. 

5.3.6. The arclike component L 

At i>lsr = —3.2 km s _1 in Fig. [HI a prominent absorp- 
tion feature is detected at approximately the location of 
the Bright Bar. This feature, indicated by red contours 
in Fig. [TU1 was included by vdWG89 in component C, 



based on the agreement in velocity. The higher resolu- 
tion provided by the present data however shows that 
this feature is distinct from component C, as will be dis- 
cussed in Sect. 17.21 This feature, which we denote by 
L, has an arc-like structure, open towards the southeast, 
similar to the shape of component D (indicated by dark 
blue contours in Fig. [TO]) . It reveals a velocity gradient 
with more negative velocities towards the southeast (e.g., 
Fig. [TB"]). similar to component D. 

5.3.7. The elongated component J 

Southeast of component L, and at approximately the 
same velocity (~ —3 km s _1 ), an elongated Hi feature 
is found in the opacity images, crossing the Bright Bar 
orthogonally. A PV diagram through this feature, which 
we label J, shows that it is kinematically distinct from 
component L, and at slightly more negative velocities. 
Component J is similar to component G in appearing 
elongated, and has a velocity gradient oriented along its 
long axis. 

5.3.8. Component K 

A compact absorption feature at ~ —3.2 km s , which 
we label component K, is detected in Fig. [HI south of 
component J. 

5.3.9. Features at the velocity of OMC-1 

Several features at the velocity of the background 
molecular cloud OMC-1 are found in the present data. 
Component H was already discovered by vdWG90 and 
is detected as an extended feature in Fig. [9] at 7.2 — 
11.1 km s _1 , and in the PV diagram in Fig. [16] In ad- 
dition, a system of compact features is found directly 
southeast of the Bright Bar and most likely associated 
with it. These features can be seen in Fig. at a ve- 
locity of 9.7 km s -1 , and we denote them collectively 
as component M. The most prominent of these is lo- 
cated at the extreme eastern edge of the Bright Bar at 
RA = 5 h 35 m 25?8, Dec = -5° 23' 57". Following the 
Bright Bar towards the southwest, several similar fea- 
tures are detected. One of these components can be seen 
in the sample spectrum shown in Fig. [S] and the H I emis- 
sion PV diagram shown in Fig. [6] 

In the Dark Bay region, a single compact H I absoption 
feature is found at the velocity of OMC-1. This compo- 
nent (component I) can clearly be seen in the opacity 
images at 11.1 and 12.3 km s _1 . It was also shown in the 
PV diagram of H I emission in Fig. 

6. H i EMISSION ASSOCIATED WITH THE ORION 
NEBULA 

The H I emission images and PV diagrams presented in 
Sect. 2] show a number of separate features revealing the 
neutral environment of the Orion Nebula and the effects 
of Hn region and the ONC on this environment. 

1. Hi emission at approximately the velocity of the 
main absorbing components A and B is detected 
in regions where strong absorption is absent. This 
layer, which is shown in the PV diagrams in 
Figs.|3HB] contains a bright Hi emission feature, de- 
tected in Fig. |4] at a spatial offset of appoximately 
150" This feature is therefore located between M42 
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Fig. 18. — Spectrum of Hi brightness temperature Tb, averaged 
over a 45" X 36" rectangular region parallel to the Bright Bar (po- 
sition angle 55°), centered on RA=5 h 35 m 21?3, Dec=-5° 25' 08". 
This region is indicated by a yellow rectangle in Fig. \3\ 

and M43 and may represent photodissociated gas 
outside an IF bounding M42 on the side of the 
Northeast Dark Lane. This reg ion lies immedi- 
ately n ortheast of sample 5-east in 10 'Dell fc Harris! 
(pOll . Their Fi gure 1 shows that this region corre- 
sponds to an overlap of a northern profusion from 
the Dark Bay and the Northeast Dark Lane. Since 
Hi emission could arise from both features, it is 
impossible to unambiguously assign the observed 
emission. 

2. A prominent Hi emission feature is detected in 
Fig. [U southeast of the Bright Bar, i.e., the Orion 
Bar PDR. As shown in Figs. [5] and [SJ an elongated 
Hi emission system extends from this feature to- 
wards the southwest at ulsr ~ 14 km s" 1 along. 

3. High- velocity Hi emission arising from a small re- 
gion is seen at ulsr ~ 18 km s _1 in Fig. [5] and in 
the PV diagrams in Figs. [6] and [Jj 

4. Hi emission is detected from the Dark Bay region, 
and this feature is centered on the H I absorption 
component I, as shown in Fig. [7] 

6.1. Hi emission from the Orion Bar PDR 

The brightest H I emission in Fig. [2] is found directly 
southeast of the Bright Bar, thus arising in the promi- 
nent edge-on PDR. An H I emission spectrum, averaged 
over the region shown by the yellow box in Fig. [31 is 
shown in Fig. 1181 This spectrum is very similar to the 
H I spectrum of this region obtained almost 40 years ear- 
lier with the Parkes 64 m telescope shown in Fig. 6 of 
IRadhakrishnan et all l| 19721) . 

Our observation of Hi emission from the Orion 
Bar probes a significantly younger PDR, with much 
higher gas density than previous studies of PDRs as- 
sociated with evolved Hli regions such as NGC281 
(|Roger fc PedladHMl . NGC1579 (|Dewdnev fe Roger! 
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Fig. 19. — Normalized crosscuts across the Bright Bar and the 
Orion Bar PDR in various tracers. The crosscuts cover a length 
of 75" perpendicular to the Bar at a position angle of 145° and 
have been averaged parallel to the bar over a strip with a total 
width of 45". The center of this strip passes through coordinates 
RA = 5 h 35 m 21?3, Dec = -5° 25' 08", i.e., the central position 
for the spectrum shown in Fig. 1181 The location of this strip is 
shown by a purple rectangle in Fi g. [3] and it is identical to that 
used by Van der Wcrf et al. (1996) to construct their Fig. 9. The 
abscissa shows distance (positive towards the southeast) from the 
IF, where the position has been chosen as the point of maximum 
3.3 fim PAH emission, which marks the location of the IF. The 
21 cm continuum and H I brightness temperature dat a are from the 
present paper, while the 3.3 fim PAH data are from lTTelens et al.l 
1T993T) and lBregman et all lfl994Tl . The H 2 v = 1-0 S(l) and C ls O 
2-1 data are from I Van der Werf et al.l 01)961 ). 

[1981 fl986h IC 5 146 ([Roger fc Irwinl Il98l . S187 
([Joncas et al.l Il992f ). S185 (which contai ns the well- 
studied PDR IC 63 IBlouin et al.l [l997h and S270 
(|Roger et al.l 120041) . Due to the high resolution and the 
edge-on orientation of the Orion Bar, our observations of- 
fer for the first time the opportunity to use H I to study 
the stratified structure of such a PDR. We therefore 
construct a crosscut perpendicular to the bar in various 
tracers, as shown in Fig. [JJJ]). This figure dramatically 
confirms the layered structure of the PDR, and for the 
first time observationally pinpoints the location of atomic 
hydrogen in a resolved edge-on PDR. The edge-on IF is 
marked by the peak in 3.3 /xm PAH emission, which is 
strongly excite d in the neutral UV-exposed layer direc tly 
outside the IF ([Tielens et al.lll993t iKassis et al.ll2006ft . 

The most important parameter in determining the 
thickness of a homogeneous PDR and therefore the sep- 
aration of the various trac ers in Fig. [19] is the dissocia- 
tion parameter y y 'nu (e.g.JTielens fc Holleiiba ch 1985b; 
Sternberg! Il988t iBurton et all 119901 : iDraine fc Bertoldil 
19961: iHollenbach fc TielenT ~[l999;) , where x repre- 
sents the intensity of the incident UV radiation field, 
and Tin is the number density of hydrogen nu- 
clei. The structure can be modified if the gas is 
clumpy (IBurton et all 119901 iMeixner fe Tielensl 119931: 
iTauber et al.lll994t lYoung Owl et al.ll2000l ). since in this 
case UV photons penetrate deep er into the clo ud through 
the interclump medium (e.g., Bois se 1990). For the 
Orion Bar PDR, the incident UV radiation field at 
the IF is x = 2.6 x 10 4 , where y is expressed in 
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units of the iDraind (|1978| ) interstellar radiation field of 
2.7 x 10 -3 erg s _1 cm~ 2 . This value of x was derived 
directly from the properties of th e Trapezium stars and 
their d istance from the Bright Bar (|Tielens fc Hollenbachl 
Il985a[) . Analysis of the stratified structure then yields 
a gas density % = 1 — 2 x 10 5 cm" 3 , (iTielens et al.l 
[1991 ISimon etallH997t [Van der Wiel et al.ll2009D . with 
higher d ensity (up to = 10 7 cm~ 3 ) embedded - 

clumps dTauber et all 119941: iVan der Werf et ail 119961: 



lYoung Owl et al.ll2000l: ILis fc Schilkdl2003D . We now dis- 
cuss our H I observations of the Orion Bar in the context 
of this model. 

As shown in Fig. [TH1 the peak H I brightness temper- 
ature in the Orion Bar PDR is approximately 230 K, 
and the H I spin temperature and kinetic temperature 
are therefore at least this high. Such temperatures 
are easily reached in PDRs exposed to an intense UV 
radiation field. For the values of \ an d «h 2 appli- 
cable to the Orion Bar, PDR models predict gas ki- 
netic temperature s of approxima t ely 1000 K at the UV- 
exposed surface (ILe Petit et al.l [2006: Mciicrin k et al.l 
l2007t Kaufman et al.lll999| ). decreasing to values of a few 
hundred K at larger distances (> 0.01 pc) from the IF. 
An upper limit for the kinetic temperature follows from 
the observed linewidth, using 



Adfwhm = 2W21n2 



kTi 



kin 



TORI 



(2) 



where tohi is the mass of the Hi atom. While the 
low velocity side of the emission feature in Fig. [JS] may 
be affected by absorption, the high velocity side ap- 
pears to be unaffected, and at that side the HWHM of 
the line is 2.6 km s , corresponding to a FWHM of 
5.2 km s _1 . This line width implies tha T s < 590 K. 
The turbulent linewidth in this region, as measured 
from the optically thin C O 3—2 line (jJohnstone et al.l 
2003) is 1.3 km s _1 . Allowing for this turbulent ve- 
locity gives a thermal FWHM of 5.0 km s _1 , im- 
plying Tki n = 540 km s _1 . Since the peak bright- 
ness temperature is 230 K, the implied peak optical 
depth is To = 0.6, so the Hi line emission is only 
marginally optically thin. The implied H I column den- 
sity is JV(Hi) = 3 x 10 21 cm -2 , averaged over the region 
used to construct Fig. [T5J The line of sight molecular 
hydrogen column density in the bar is iV(H2) k 6.8 x 
10 22 cm" 2 (IHogerheiide et al.l Il995t IJansen et al.l Il995t 
IVan der Werf et al.l 119961: I Van der Wiel et al.l I2009D . so 
that only about 2% of all hydrogen nuclei are in atomic 
form. Since the dense clumps in the PDR will likely re- 
main molecular, it is more appropriate to consider only 
the interclump medium with nn = 1 — 2 x 10 5 cm~ 3 
for calculating the local Hi abundance. The line-of- 
sight depth of the IF in the Bar region i s approximately 
0.12 p c, as derived from modeling by iPellegrini et al.l 
(2009). Using these numbers, the Hi abundance in the 
Bar is 5 — 10% at the position of the Hi peak. 

The derived atomic fraction indicates that the observed 
H I emission originates in the region where the transi- 
tion towards molecular gas occurs. This conclusion is 
supported by the H2 temperature in t his region from ro- 
tational lines, wh ich is 400 - 700 K (jAllers et al.l 120051: 
iShaw et aI1l2009D . in excellent agreement with our esti- 



mate of 540 K for the H I. 

Inspection of Fig. [19] reveals strong H 1 emission from 
the region where H2 vibrational line emission shows 
a maximum. This agreement is physically significant. 
Direct H2 photodissociation from the ground state is 
strongly forbidden, since the molecule is homonuclear. 
The actual photodissociation of H 2 is a two-step pro- 
cess, initiated by the absorption of UV photons in the 
Lyman or Werner bands, as first pr oposed in 1965 b y 
Solomon (private communication in iField et"aTI 119661 ) . 
The UV absorption is followed by a radiative cas- 
cade, in which there is an 11% chance of dissociation 
(jStecher fc Williams! [1961 . In the remaining 89% of 
cases the molecule cascades down to the ground state, 
and the resulting fluoresc ent photon emission produces 
the H2 vibrational lines dBlack fc Van Dishocck 119871 : 
lSternberd[T98l iSternberg fc Dalgarndll989l e.g~7). The 
strong Hi emission from the region of maximum H2 vi- 
brational line emission therefore directly supports the 
photodissociation mechanism, and our data reveal this 
agreement for the first time. 

At larger distances from the IF, the H 1 brightness tem- 
perature decreases slowly. The optically thin C ls O J = 
2—1 emission peaks at a larger distance from the IF than 
the H 1 emission (21" or 0.05 pc). The brightest H 1 emis- 
sion is thus located in the region where CO is photodisso- 
ciated. In this region the gas-phase carbon is singly ion- 
ized, and detected through the fC III 158 /im fine s truc- 
ture line (jStacev et al.l ll993: Herrm ann et al.lll997f) and 
recombination lines in the radio ( Jaffe fc Pankoninlll978l: 
INatta et al.lll994HWvrowski et al.lll997D and in the near- 
infrared ([Walmslev et al.l l2000f) . The slower decrease of 
the Hi brightness temperature towards larger distances 
from the IF compared to H2 v = 1—0 S(l) results from 
the fact that co llisional excitation contrib utes to the flux 
of this H2 line ([Van der Werf et al.| [l996'). Observations 
of the fainter H2 v = 2—1 S(l) line, which is dominated 
by UV- p umped fluorescence r eveal a slower decline from 
the IF (lHavashi et al.l 119851: IVan der Werf etldl 1199 ' 



Luh man et al.lll998t iMarconi et al.lll998t IWalmslev et al 
2000), matching the decreasing H 1 brightness tempera- 



ture in the same region. 

The Orion Bar arises from an escarpment protruding 
from the OMC-1 molecular cloud towards the observer. 
Southeast of the Br ight Bar the IF cur ves back to a 
more face-on aspect (|Wen fc Q'Delll [l995). This geome- 
try is illustrate d schem atically for inst ance in Fig. 3 of 
IPellegrini etail ([200l and Fig. 13 of IQ'Dell fc Hamsl 
d2010(). The~present H 1 emission results show that the 
northeast-southwest extent of this escarpment is much 
larger than the length of the bright IF in the Huygens 
region. As the highly elongated Hi emission feature at 
u lsr = 14-9 km s~ 1 (Fig. [2]) shows, the escarpment ex- 
tends significantly beyond the Bright Bar towards the 
southwest. This Hi emission feature can be seen as a 
separate velocity component at ulsr ~ 15 km s _1 in 
Fig. [SI extending over about 11' (1.5 pc). This Hi feature 
is located in the EON, which is much less well studied 
than the bright Huygens region. The extended Hi fea- 
ture has a counterpart in a similarly extended feature ob- 
served in 8 /im emission using IRAC on the Spitzer Space 
Telescope, representing emission from PAHs at 7.6, 7.8 
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and 8.6 ixirj^. An IF bounding the elongated H I feature 
on its northern side is detected in the wide-field opti- 
cal images of the Or ion Nebula obtained with ACS/HST 
(jHennev et al.l [2007b . This geometry confirms that the 
Orion Bar PDR, and thus the escarpment from OMC- 
1, extends significantly beyond the bright section in the 
Huygens Region. 

Southeast of the brightest section of the Orion 
Bar P DR the IF cu r ves back to a more face-on 
aspect (IWen fe O'Delll [l995t IHoprheiide et all [19951: 
Jansen et al.1119 95; Pellegrin i et al.H2009t lAscasib ar et al.l 
20111) . A recent Spitzer IRS study has revealed [Nelil] 
15.56 fim and other ionized gas lines out to 12' south- 
east of the Trapezium stars, i.e., far beyond the Bright 
Bar ([Rubin et al.ll20H"l ). Our data reveal an extended 
Hi cloud in this region with a central velocity vlsr ~ 
13 km s . This feature may represent atomic gas in the 
extended PDR with a more face-on aspect here than in 
the bright Bar region. All details of the geometry of the 
Orion Bar region are represent ed in the diagram show n 
in the upper panel of Fig. 13 of lO'Dell fe Harris! (|2010l ). 

6.2. Hi emission from the Dark Bay region 

The H I emission images in Fig.[5]also reveal a ~ 45" di- 
ameter Hi emission feature at velocities up to 17 km s _1 
located in the Dark Bay area. The peak of this fea- 
ture is at RA = 5 h 35 m 24?0, Dec = -5° 22' 30" at 
14.9 km s — the feature extends towards the southeast 
where it crosses the compact Hi absorption feature I at 
RA = 5 h 35 m 25?4, Dec = -5° 22' 38". Given the excel- 
lent agreement both in velocity and position, as shown 
in Fig. the emission and absorption feature are almost 
certainly physically related. 

Given that the continuum brightness temperature of 
the H II region in this area is T c ~ 1500 K, the emitting 
H I must be located behind the ionized gas, and is there- 
fore not associated with the Dark Bay, which represents 
a tongue of absorption in front of the Hn region, with 
a velocity i>lsr ~ 6 km s _1 as measu red from radio re- 
comb ination lines of partly ionized gas ( Jaffe fc Pankoninl 
[19781) . 

The H I emission in the Dark Bay area is located east of 
the region where th e IF curves to a more edge-on orien- 
tation, as shown by IWen fc O'Delll (|1995l) . This location 
suggests a geometry analogous to that of the Bright Bar, 
discussed in Sect. 16.11 This idea is supported by the 
presence of a 13 CO J = 3—2 emission feature in this re- 
gion, approximately coinciding in or ientation and exten t 
with the Hi emission (see Fig. 12 of lBuckle et al.|[2012T) . 
Dust emission from this regi on has been detected wit h 
SCUBA at 850 and 450 /xm (| Johnstone fc Ballvl[T999l) . 
Thus the H I emission in the Dark Bay may originate in 
an approximately edge-on PDR located behind the Hn 
region. The H I absorption feature I, which is associated 
with the emission, most likely represents a dense region of 
this PDR. The expanding IF will propagate more slowly 
where it encounters dense neutral gas, creating a concave 
region in the approximately edge-on IF. The H I in this 
concave region then has a background radio continuum, 
giving rise to a detection in absorption. The dust in com- 
ponent I will obscure the section of the IF lying behind 

8 http:/ /www. spitzer. caltech.edu/images/1648-ssc2006-16b- 
The-Sword-of-Orion 



it (as seen from Earth). Therefore this model also pre- 
dicts enhanced extinction at the position of the absorbing 
clump. As shown in Fig. 1101 this enhanced extinction is 
in fact observed, since Hi component I corresponds accu- 
rately in pos ition, extent and orientation w ith extinction 
"Knot 2" of lO'Dell fc Yusef-Zadehl (pOOOt) . 

6.3. High velocity Hi emission from the 
Kleinmann-Low region 

High positive velocity H I from a region approximately 
45" in diameter (e.g., at ulsr = 18-7 km s _1 in Fig. [2]) 
extends to LSR velocities of approximately 31 km s _1 , 
as shown in Figs. [5] and [7] An image of this Hi emis- 
sion, integrated over the LSR velocity range from 18.7 
to 31.0 km s _1 is shown in Fig. [20l The high velocity 
H I emission coincides in position with the high velocity 
outflow in the Orion-KL region, embedded in the OMC-1 
cloud behind the optical Huygens region. 

The high velocity outflow in the Orion-KL region ex- 
hibits a wide opening angle and line- of-sight velocities 
of up to 100 km s _1 in CO lines (jZuckerman et all 
[19761 IKwan fc ScoviUd 119761) . The shocked gas in this 
system produces luminous H2 vibrational line emission 
(|Gautier et al.l 119761: INadeau fc Geballd 119791 ). with a 
spectacular morphology displaying multiple "bullets" 
(bow-shoc k tips) and "fingers" ( bow-shock wakes) a s 
shown bv lAllen fc Burtonl (|1993D . iStolovv et all (|1998l) . 
Schul tz et al.l ([1999D . and illustrated in the H 2 v = 1— 
S(l) image shown in Fig. [501 High resolution CO ob- 
servati ons reveal similar high ve l ocity "bullets" and "fin- 
gers" Jphernin fc Wrightl 119961 : iRodriguez-Franco et al.l 
1999bJaT lZapata et al.ll2009D . emerging from a common 
position. A number of optically ide ntified HH ob- 
jects originate from the same position (|Doi et al.l 120021 : 
lO'Dell fe Hennevll2008l : IZapata et aTll2009l ). 

The lack of correspondence of the high velocity Hi 
emission with prominent H2 and CO "fingers" , and the 
absence of Hi emission at velocities i>lsr > 31 km s , 
indicates that the Hi is not associated with the high- 
est velocity outflowing gas. Nevertheless, the velocity of 
the observed Hi with respect to OMC-1 shows that the 
material does participate in the outflow. The two high 
velocity features separated by a local minimum (Figs. [S] 
and [7J suggest that both outflow lobes are detected, al- 
though the northwest lobe is much more prominent in H I 
emission than the southeast lobe (Fig. |20|) . The north- 
west lobe displays both redshifted and blueshifted Hi 
emission, and this situation matches that observed i n 
CO (e.g., IChernin fc Wrighd fl99i IZapata"et~aT] [20091) . 
The southeast lobe is much less well-defined (both in H 1 
and in CO), and it has been suggested that the outflow 
towar ds the southeast is blocked by the Orion-K L Hot 
Core (jChernin fc Wrightl fl99fit IZapata et al.ll2011bh . 

Figure [20] shows a positional match between the high 
velocity H 1 emission and the shocked H2 vibrational line 
emission, but no detailed correspondence. In fact a de- 
tailed correspondence is not expected since Fig. [20l shows 
only the redshifted Hi, which is moving into OMC-1, 
away from the observer. The H2 emission on the other 
hand shows all shocked H2, regardless of velocity, but 
may be biased in favor of blueshifted gas with lower ex- 
tinction. The H 1 is clearly concentrated in the region of 
the flow closest to the principal heating sources of the 
Orion-KL region. The lack of more extended Hi emis- 
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Fig. 20. — Left panel: image of Hi emission associated with Orion-KL, integrated over 18.7 < Ulsr < 31- km s -1 (grayscale with white 
contours). Black contours indicate 1420 MHz continuum surface brightness levels of to 0.5 mjy be am -1 in steps of 0.0 4 mjy beam -1 . 
The cyan circle indicates the origin of the high velocity outflow as derived from the CO observations {Zapata ct al. 2009). The dark blue 
stars represent the Trapezium stars. The resolution of the Hi data is indicated by the 7'.'5 circular beam in the lower righthand corner 
of the diagram. R ight panel: the same H I emission image represented by black contours, superposed on the H2 v = 1—0 S(l) i mage of 
I Bally e t al. (2011). The cyan dot is as in the left panel, and the blue dot represent the position, taken from Goddi ct al. (2011), of the 
Becklin-Neugebauer object. 



sion may result from a decreasing column density as the 
flow expands away from its source, but may also repre- 
sent photodissociation of the molecular material by UV 
radiation close to the heating sources. 

The high velocity Hi matches well, both in size and 
velocity, with the e xpanding CO shell re cently discov- 
ered in Orion-KL bv lZapata et al.l (|2011a[ ). and centered 
approximately on the outflow origin. This bubble has a 
diameter of about 25", an expansion velocity of about 
10 km s -1 , and a dynamical age of 500 — 1000 years. 
The left panel of Fig. [20] indeed suggests that the ori- 
gin of the outflow is located at a local column density 
minimum, bounded north, east and west by an incom- 
plete and clumpy shell. Inspect ion of the CO PV diagram 
shown by Zap ata et al.l (|2011aD shows that the northwest 
region of the CO shell is prominent in both redshifted and 
blueshifted emission, while in the southeast part the red- 
shifted emission dominates, in agreement with the sit- 
uation observed in Hi. The precise relation of the ex- 
pending CO shell to the rest of the outflow system is 
unknown. However, our high velocity Hi data show fea- 
tures of both the CO shell and the inner region of the H2 
emission, suggesting that these are directly related, and 
plausibly trace the same gas. 

The mass in redshifted Hi represented in Fig. [20] 
is only 2 x 10 -5 M Q . Based on the limited velocity 
range sampled, the total mass of Hi associated with 
the outflow could well be a factor of 3— 5 higher, i.e., 
(0.6 - 1.0) x 10~ 4 M© (0.06 - 0.1 Jupiter masses). This 
mass is insignificant compared to the total mass of molec- 
ular gas i n the expanding bub ble (~ 5 M©) as derived 
from CO (|Zapata et al.ll2011al ). 



7. H 1 ABSORPTION ASSOCIATED WITH THE 
ORION NEBULA 

The H 1 absorption dataset presented in Sect. [5] shows a 
striking number of structures. Nevertheless, the various 
velocity features can be subdivided into a small number 
of groups, according to common properties: 

1. large-scale features covering M42 as well as M43. 
These are the components A and B, discussed by 
vdWG89. 

2. arclike features. These are components D and F of 
vdWG90 and L identified in the present dataset. 
in addition, component C of vdWG89 can be rec- 
ognized as an extended arc; 

3. elongated features displaying unique kinematic sig- 
natures. In the present data, these are components 
G and J. 

4. features at the velocity of the background molecular 
cloud. These components are H, I and M. 

Based on the agreement in velocity with the ionized 
gas, which is streaming away from OMC-1 towards us, 
vdWG89 interpreted component C as neutral material 
closely associated with the ionized gas, and suggested 
that this material was entrained in the ionized flow. For 
the blueshifted small-scale components, vdWG90 pro- 
posed acceleration by the rocket effect (jOort fc S pitzcr 
UMB) as the origin of the blueshift with respect to the 
bulk of the neutral gas. For these explanations to be 
valid, ionization fronts would have to be present on the 
surfaces of the Hi features. However, HST and ground- 
based observations have failed to reveal optical rims 
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associated with these compo nents ([Hennev et al.l [20071 
IGarcia-Diaz fc Hennevl[2007f) . Therefore, new interpre- 
tations for these features will be developed in the follow- 
ing subsections. 

7.1. The large-scale atomic Veil 

The large-scale structure of the atomic Veil covering 
M42 and M43 is dominated by Hi velocity components 
A and B, and the red wing of the absorption lines is de- 
tected at 7.2 km s _1 < u LSR < 9.7 km s _1 in Fig. El 
In agreement with vdWG89, a strong opacity gradient 
found towards the Northeast Dark Lane; a positive ve- 
locity gradient of about 2.4 km s _1 pc _1 is found in 
the same direction. This velocity gradient is detected 
in both components, indicating that they are related 
(vdWG89). vdWG89 suggested that component A rep- 
resents Hi in the envelope of OMC-1, extending in front 
of M42 and M43. A small velocity difference (H I compo- 
nent A is blueshifted with respect to the molecular gas 
by ~ 2 km s _1 ) represents a slow expansion of the cloud 
envelope. 

7.2. Arc-like Hi features 
7.2.1. The southwest arc: component C 

The Hi absorption feature at «lsr ~ —3.2 km s _1 , 
forming an almost semicircular arc in the southwest part 
of M42, is Hi component C of vdWG89. The opacity 
images at ulsr = —4.4, —3.2 and —1.9 km s _1 show 
that this feature contains a string of clumps. 

The system of optical flows in the southern part of 
M42 originates mos tly from a reg i on on the east side of 
the Orion-S region; iHennev et al.l (|2007f ) referred to this 
location as the Optical Outflow Source (OOS; indicated 
as such in Fig. [T0| . Molecular outflows o riginate from 
the O rion-S cloud (summarized in Fig. 1 of lO'Dell et al.l 
2009). Flows originating in this area branch out in mul- 
tiple directions, driving several well-k nown HH objects 
([Bally et all 120001 : IHennev et al.ll2007h . The arc formed 
by Hi component C opens towards Orion-S. This feature 
may thus result from a flow originating from the Orion-S 
region. 

A flow system located at the position angle from Orion- 
S towards component C was detected as an extended 
arc of optical em i ssion, denoted as the Southwest Shock 
bv IHennev et al.l (|2007l . see their Fig. 6). This feature 
is located approximately twice as far from the OOS as 
the arc formed by component C, but its shape is very 
similar to the Hi. These features may therefore result 
from two separate ejections from the same object in the 
OOS. However, further data would be needed to verify 
this hypothesis. 

The PV diagram shown in Fig. [12] shows the relation 
between component C and components A and B, that 
define the large-scale structure of the Veil. At small neg- 
ative spatial offsets (i.e., towards the northeast), both 
components A and B are detected, but component C is 
not observed. For instance, at offset ~ —30", where the 
lines are not saturated, both components A and B are 
clearly detected as kinematically separate components. 
This situation is different at positive offsets. Where com- 
ponent C appears (at offset ~ 15"), component B disap- 
pears. This result clearly suggests that the arc formed by 
component C consists of material swept up from compo- 



nent B, displaced towards a more negative velocity. The 
fact that component B is not detected in the region inside 
the curved arc (Fig. 1111) is consistent with this interpre- 
tation. 

7.2.2. The arc-like component D 

Component D was first identified by vdWG90 and 
is the largest of the small-scale components, extend- 
ing over a significant part of the Huygens region be- 
tween the Trapezium stars and the Bright Bar. As 
noted by vdWG90, this component coincides in posi- 
tion and velocity with the region where velocity split- 
ting is observed in several optical emission lines from 
M42. This velocity splitting was first studied in detail 
bv iDeharvengl (|1973f) . whose [Nil] 6583 A line-splitting 
region A corresponds to our Hi component D. A full 
ki nematic atlas of several optical lines has been presented 
bv IGarcia-Diaz fc Henn ev (2007) and Garcia-Dia z et al.1 
(2008) ; this dataset shows the line splitting in detail. The 
ionized gas component matching our H I component D 
is referred to as the "Southea st Diffuse Blue Layer" by 
IGarcia-Diaz fc Henneyl ([20071 ). who note that this com- 
ponent is detected in [Sn] 6716 and 6731 A but not in 
[S III] 6312 A. This result is important since it implies 
that the ionizing spectrum is rather soft , and not pro- 
vided by 9 1 C Ori. IGarcia-Diaz fc Hennevl (|200l suggest 
that the blueshifted velocity component is ionized by the 
star 9 2 A Ori, located southeast of the Bright Bar. 

The present results shed new light on this issue. The 
arc-like morphology of component D was already pointed 
out in Sect. 15.3.21 The kinematic structure of this com- 
ponent suggests an expanding shell with center of ex- 
pansion in the direction towards which the arc opens, 
i.e., southeast of the Bright Bar. Figure [T0l shows that 
the arc opens towards the star 9 2 B Ori, located approx- 
imately on it s axis of symmetr y; this star has spectral 
type B0.7V (|Sim6n-Diazl 12010ft . The observed geome- 
try suggests that this star provides the ionization for the 
Southeast Diffuse Blue Layer. In this model Hi com- 
ponent D represents the neutral material outside of this 
Hn region. This explanation implies that an IF should 
be present between the Southeast Diffuse Blue Layer and 
H I component D; the presence of an IF is confirmed by 
the detection o f weak [O i] 6300 A em ission associated 
with this layer ([Garcia-Diaz et al"1l2008[ ). 

7.2.3. The arc-like component L 

The Hi component L at i>lsr = —3.2 km s _1 (shown 
by the red contours northwest of the Bright Bar in 
Fig. [TO)) displays an arc-like structure similar to that of 
component D, also opening towards the southeast. The 
characteristics of expansion are seen in Fig. I13| where it 
exhibits a velocity gradient with more negative veloci- 
ties towards the southeast. For this component, the star 
9 2 A Ori is located close to the symme try axis of the arc. 
The s pectral type of this star is 09V ([Simon-Diaz et all 
2006). This configuration indicates a model identical to 
that described above (Sect. I7.2.2|) for component D, ex- 
cept that here 6 2 A Ori is the exciting star. 

7.2.4. The expanding shell component F 

The third component to display an expanding shell-like 
morphology is component F. This component is located 
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close to a large complex of extended bow shocks in the 
ionized gas related to the Herbig-haro object HH202. 
As can be seen in Fig. [TOl HH 202 is located (in pro- 
jection) at the western edge of the arc defined by Hi 
component F. At this position, component F splits into 
two separate velocity components (already recognized by 
vdWG90 and indicated Fl and F2 in that paper). The 
velocity splitting is evident in Fig. 1151 where compo- 
nent F is observed to split into two subcomponents with 
different spatial and kinematic offsets. The subcompo- 
nent at the most positive spatial offset matches well in 
position with HH 202. 

The positions of high proper motion features in the 
ionized gas, shaped liked parabolic a rcs, have been es- 
tablished by IQ'Dell fc Hermevl (120081) . using HST im- 
ages from several epochs. The driving source for HH 202 
is located in the O OS region associated with Orion-S 
(jHennev et al.ll2007[ ). Indeed, a jet connecting HH202 to 
this region has been detected in th e shock-tracing [Fen] 
1.257 (im line (jTakami et al.ll2002D . 

Two explanations can be considered for the presence 
of atomic hydrogen associated with, but not precisely 
superposed on HH 202. 

1. It is possible that HH 202 and H I component F are 
formed by a jet creating a bow shock in both the 
ionized gas (observed as HH 202) and the neutral 
Veil (obs erved as Hi compo nent F), as proposed 
earlier bv IQ'Dell et all ffT997h . 

2. Alternatively, it is possible that Hi component F 
results from rapid recombination in the dense post- 
shock gas associated with HH202. This mode l 
has been proposed by Mcsa-Dclgado et al. (2009), 
who show that in HH202— S (the brightest knot in 
HH 202) the post-shock density is sufficiently high 
to trap the IF. As a result, the gas behind this 
dense post-shock region is shielded from the ion- 
izing radition of the Trapezium stars, and rapidly 
recombines and cools. 

3. Finally, it is possible for this component that the 
acceleratio n is provided by the rockets effect de- 
scribed bv lOort & Spitzerl (|1955f ). 

In the second and third models, the neutral Hi should 
be located behind the IF in the direction away from the 
Trapezium stars, i.e., west of the bright rims of HH202. 
However, none of the absorbing H I is located in this re- 
gion, and most of it is in fact located east of HH202. 
This geometry argues against the last two models. In ad- 
dition, a strong argument in favour of the first of these 
three models comes from the fact that the blueshifted 
H I gas forming component F corresponds to a gap in the 
Veil, as observed in Fig. [15] This is exactly the geometry 
that is expected when part of the Veil gas is accelerated 
to a more negative velocity by the impact of a jet. In 
summary, the location of the H I with respect to HH 202, 
and its detailed kinematic structure argue in favour of 
an interaction of the jet driving HH 202 with the neutral 
Veil. 

7.3. Elongated Hi features 

Among the kinematic H I features identified, compo- 
nents G and J have an obvious elongated appearance, 
quite distinct from the arc-like features discussed above. 
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Fig. 21. — Hi opacity contour s overlayed on the H2 v = 1—0 
S(l) image of Bally et al. (2011). Light blue contours show Hi 
opacity at ulsr = —9-6 km s _1 from Hi components Gl and 
G2, with contour levels from 0.1 to 0.6 in steps of 0.1, then from 
0.8 to 1.8 in steps of 0.2. Yellow contours show Hi opacity at 
^LSR = — 0.6 km s — 1 from component G3, with contour levels 
from 0.50 to 2.25 in steps of 0.25, then from 2.5 to 4.0 in steps of 
0.5. 

7.3.1. Component G 

The Northern part of M42, where component G is 
located, harbors a system of optical HH objects (e.g., 
HH201, HH 205-210, all indicated in Fig. [10]), which 
are associated with the complex and spectacular sys- 
tem of "fingers" observed in H2 vibrational line emis- 
sion (lAllen fc Burtonlll993HSalas et al.lll999l : lB~ariv et al.l 
120111 ) and CO (jZapata et all l2009h . Emission in the 
shock tracing [Fen] 1.64 /xm line is found at the tips 
of the fingers. The fact that these "fingers" display op- 
tical line emission shows that in this region they have 
bro ken out of the bac kground molecular cloud towards 
us (|Q'Dell et al.l[l997l) . 

In the present H 1 data, several velocity systems are 
found in this region, as can be seen in Fig. 1171 The 
highest (negative) velocities are found for component G2 
(light blue contours in Fig. [TU|) . and this component is, 
in location as well as orientation, closely aligned with 
the line connecting the "dynamical center" where this 
outflow system originates (as discussed in Sect. I6.3|) to 
the Herbig-Haro objects HH 205—207. This line also cor- 
responds to one of the most prominent molecular "fin- 
gers": the 1 : 00 sy stem in the (hour dial) notation of 
IQ'Dell et ail (|1997ft . This agreement is shown in Fig. EH 
where component G2 is shown in the blue contours at 
Declination north of -5° 21' 40". South of this Declina- 
tion, the blue contours trace component Gl, and Fig. |2~T1 
shows that this feature does not correspond to any of the 
shocked H2 features and, given its orientation, does not 
originate in the Orion-KL region. At less negative ve- 
locities, component G3, which is approximately aligned 
with Gl but displaced eastwards, is prominent, as shown 
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in Fig. HUby the yellow contours. If components Gl and 
G3 originate from a common source, this source must be 
located approximately 30" east of the Trapezium stars. 
However, this region does not harbour any known driving 
source for these flows. 

The asssociation of component G2 with the H2 "finger" 
system related to the Orion-KL outflow enables a new es- 
timate of the line-of-sight location of the dynamical cen- 
ter of the outflow behind the IF. The Herbig-Haro object 
HH 205 is associated with this H2 "finger" , and closely 
aligned with H 1 component G2 (Fig. ITO)) . The radial ve- 
locity of HH205 is v LS -R = -200 km s^ 1 (|Q'Dell et al.l 
Il997t ). while its ta ngential velocity has been measured 
to be 320 km s" 1 (jBallv et al.ll2000f) . This system thus 
follows a trajectory with an angle of ~ 30° relative to 
the plane of the sky. The position where the H2 "finger" 
breaks through the IF is displaced from the dynamical 
center of the outflow by a projected distance of 0.14 pc. 
Combining these results yields a distance of the dynami- 
cal center of the outflow of ~ 0.1 pc behind the IF. This 
dis tance is somewh at smaller than the ~ 0.2 pc derived 
bv lDoi et all (pOOl using a similar calculation based on 
the velocity vector and location of the Herbig-Haro ob- 
ject HH201 (jGraham et al.ll2003h . 

7.3.2. Components J and K 

The other elongated H 1 feature in our dataset is com- 
ponent J, which crosses the Bright Bar orthogonally. 
This feature is best seen at klsr = —3.2 km s _1 in 
Fig. [5J Both in position and orientation, this feature 
closely coincides with the large, bright pair HH 203/204. 
Inspection of the H 1 opacity images in Fig. [9] shows a 
peak in H 1 opacity located at the position of HH 203 at 
LSR velocities from —7.0 to —4.3 km s . However, at 
—3.2 km s _1 an extension towards HH204 is apparent, 
so that Hi is assoc iated with both Herbig-Haro objects. 
lO'Dell et al.l (|1997l ) suggested that these objects mark 
the points where collimated flows impact the extended 
H 1 Veil. The detection of H 1 associated with these fea- 
tures provides support for this model. 

Hi component K, located south of J, coincides 
with a region rese mbling a bow shock (see Fig. 6 of 
iHennev et all 120071 ) . This region displays high ioniza- 
tion emission ([Oni]) as well as prominent lower ioniza- 
tion ([Nil] and [Sil]) ridges. This region may represent 
an outer shock related to the HH 203/204 system and 
would again be the result of a collimated flow striking 
the Veil. 

7.4. The blueshifted H I feature E 

The remaining blueshifted H 1 feature, component E, 
appears neither arc-like nor elongated and thus does not 
fit into any of the above categories. 

Component E is located in the region of the Dark Bay. 
This feature consists of 10 unresolved or barely resolved 
clumps, forming a coherent velocity structure fFig. [14]). 
Due to the lack of detailed optical information in the 
Dark Bay region, no associations with ionized gas flows 
can be established. However, the different morphology 
of component E as well as the larger linewidth com- 
pared to the other H 1 components suggests that compo- 
nent E does not arise from the same process. As noted 
in Sect. 16.21 the Dark Bay gas is at vlsr ~ 6 km s — x . 



Component E is therefore signficantly blueshifted with 
respect to the bulk of the Dark Bay gas. An attractive 
hypothesis is that Hi component E arises in a PDR re- 
lated to the Dark Bay, at the side facing the ONC. The 
thermal expansion of the H11 region may then account 
for the blueshift of H 1 absorption with respect to rest of 
the Dark Bay gas. 

7.5. Hi features at the velocity of the background 
molecular cloud 

Hi components H, I and M have central velocities 
wlsr = 9 — 10 km s , corresponding to the velocity 
of OMC-1. Since these features are observed in absorp- 
tion, a significant radio continuum must arise behind 
them. Component I has already been discussed above 
(Sect. 16. 2[) . in connection with the Hi emission associ- 
ated with that component. 

7.5.1. Component H: Hi absorption from Orion-S 

The H 1 absorption component H (which was already 
noted by vdWG90) corresponds closely in location, ve- 
locity, extent and morpholo gy with an H2CO ab sorp- 
tion feature first identified bv lJohnston et all (11983f) an d 
studied at higher resolution by iMangum et al.l (p.993), 
as indicated in Fig. [TO] This feature is related to the 
Orion-S molecular core. The detection of Orion-S in 
absorption indicates that the IF separating M42 from 
OMC-1 must be located behind Orion-S. The implica- 
tion is that the Orion-S m olecular core must b e located 
within the ionized nebula. 10 'Dell et al.l (|2009D were the 
first to argue for this geometry, based on the detection 
of H2CO absorption and anomalies in the derived ex- 
tinction; the present H 1 data confirm this picture. The 
detailed 3-dimensional structure and the physical associ- 
ation of th e various features in this region are discussed 
in detail by 10 'Dell et al.l (|2009f) . These authors also con- 
struct a physical model of the region, showing that the 
density and column density are sufficient for Orion-S to 
survive in the harsh radiation environment of M42. 

7.5.2. Hi absorption associated with the Bright Bar 

H 1 produced in the Orion Bar PDR will be at i>lsr = 
10 km s _1 , as d erived from molecular lines in this region 
(e.g., Fig. 8 of iVan der Werf et al.l fl99l . The collec- 
tion of clumps which we have collectively labeled compo- 
nent M (as indicated in the opacity image at 11.1 km s _1 
in Fig. [9]), closely trace the outside of the Bright Bar and 
most likely represent clumps of photodissociatcd H 1 in 
the Orion Bar PDR. The small but non-zero angle of the 
IF with respect to edge-on causes most of the photodis- 
sociated gas to be located behind the radio continuum 
of the Bar; this Hi cannot be dectected in absorption. 
However, since the IF will not be flat, a background ra- 
dio continuum will be available over a small distance to- 
wards the southeast in regions where small amounts of 
H 1 are located in concave sections of the IF. Towards 
this continuum the H I can be observed in absorption. 
At larger distances from the IF, photodissociated Hi in 
the Orion Bar PDR is observed in emission, as discussed 
in Sect. O 

7.6. Non-detection of Hi from proplyds 
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No Hi absorption was found associated with known 
proplyds in the Orion Nebula, with upper limits of 
about a Jupiter mass. This result indicates that evap- 
oration and a blation of molecular gas associated with 
the proplyds (IChen et al.l [l99l IHennev fc OT)el fl999l: 
IHennev et al.l I2002T ) is rapidly followed by photoioniza- 
tion. Alternatively, the external IF confining the pro- 
plyds may be opaque in the 21 cm continuum, prohibit- 
ing the detection of neutral hydrogen. 

8. THE VEIL OF ORION 
8.1. Location of the Veil 

Several of the blueshifted velocity components show 
clear evidence of interaction with the Veil. This behavior 
is most obvious for component F, where the Veil shows a 
gap corresponding closely in position with the blueshifted 
gas of component F (Sect. I5.3.4j) . Similar evidence is 
found for component fSect. 1772.11) . For component G 
the evidence for interaction with the Veil is weaker, but 
it clearly reveals interaction with a layer of neutral gas. 

Component F therefore plays a crucial role in analy- 
sis of the interaction of the ionized flows with the Veil. 
This component is driven by the prominent HH object 
HH202, for which the space velocity vector is known, 
since both r adial and tang e ntial velocities have been 
determined (iDoi et al.l l2002l 12004 IHennev et all 120071: 
IQ'Dell fc Hennevl 120081) . This flow system originates in 
the OOS associated with Orion-S; therefore its radial dis- 
tance with respect to Orion-S can be calcula ted. This 
distance is found to be 0.17 pc (|Doi et al.ll2004fl . The lat- 
eral extent of the HH 202-component F system is about 
50" or 0.11 pc, and this size can be used as the maximum 
distance between HH202 and Hi component F. We con- 
clude that Hi component F, and therefore the neutral 
Veil, is located at most 0.3 pc in front of Orion S. 

This distance is to be compared to the ~ 1 pc dis- 
tance bet ween the Trapeziu m stars and the Veil, esti- 
mated by lAbel et al.l (|2004f ). The relative line-of-sight 
distance s of Orion-S and t he Trapezium stars are un- 
known. IQ'Dell et alj (2009) argue that Orion-S cannot 
be much behind the Trapezium stars, or else it would 
not be possible to distinguish it from main IF separating 
M42 from OMC-1. On the other hand, if Orion-S were 
significantly in front of the Trapezium stars, it would be 
in front of almost all optical line emission, and create a 
conspicuous extinction feature. While an extinction fea- 
ture at this pos ition has indeed been detect ed, it is not 
very prominent (jO'Dell fc Yusef- Zadch 2000h , and signif- 
icantly less pronounced than for instance the Dark Bay, 
which is indeed located in front of the ionized nebula. 
Considering the large column density of Orion-S, its rel- 
atively modest derived extinction must be due to the fact 
that a significant fraction of the nebular emission arises 
in front of it. Therefore Orion-S is unlikely to be much 
closer to the Veil than the Trapezium stars. 

An additional argument can be found in the morphol- 
ogy of Hi component H, which corresponds to Orion-S. 
The absorbing H I displays an opacity peak at its north 
edge, i.e., approximately facing the Trapezium stars. At 
the same position a small local maximum in radio contin- 
uum is found. This geometry suggests that the ionizing 
and dissociating radiation has created an IF and PDR 
at the side of Orion-S facing the Trapezium. The impli- 



cation is that this radiation must come mostly from the 
side, and that the line of sight connecting Orion-S and 
the Trapezium stars is unlikely to be inclined more than 
45° from the plane of the sky. This argument suggests 
that Orion-S is unlikely to lie more than 0.1 pc in front 
of the Trapezium stars. The surface of the Veil facing 
the ionized nebula will then be at most 0.4 pc in front of 
the Trapezium stars. 

Observations of UV absorption lines reveal local den- 
sities in the Veil, also providing limit s on the distance 
of th e Veil from the Trapezium stars (jAbel et al.l 120041 . 
l2006h . An upper distance limit of 3 pc was determined 
by the absence of strong H2 UV absorption lines, which 
implies a significant UV field to keep the molecular hy- 
drogen fraction at the observed value. A lower limit to 
the dista nce of the Veil to the Trapezium stars was pro- 
vided bv lAbel et alj (|2004[ ). These authors show the de- 
tailed dependence of a large number of observables on 
the Veil- Trapezium distance. In addition they used lim- 
its on H,3 and [Nil] 6583 A emission from the surface of 
the Veil facing the Trapezium to derive a minimum Veil- 
Trapezium distance of 0.33 pc (for H/3) and 1 pc (for 
[Nil]). The [ Nil] emission was however considered in 
more detail by I Abel et alj ([2006), who demonstrated the 
presence of a density-bounded layer of ionized gas located 
between the Trapezium and the Veil and moving towards 
the Veil. These authors conclude that this layer absorbs 
most of the ionizing radiation from the Trapezium stars, 
which accounts for the absence of a pronounced IF at 
the surface of the Veil facing the Trapezium. This re- 
sult has the important implication that the emission line 
surface brightness of the illuminated side of the Veil can- 
not be used to constrain its distance to the Trapezium. 
In summary, we conclude that the distance between the 
Veil and the Trapezium of 0.4 p c derived in this paper 
is consistent with the analysis by I Abel et alj (|2004f ) and 
lAbel et all ((2006T) . 

8.2. Geometry of the Orion Nebula and the Veil 

Our results reveal direct interaction of ionizing flows 
originating close to Orion-S with the Veil. The evidence 
is particularly strong for the HH202 (Hi component F, 
Sect. 17.274")) , which shows that only velocity component B 
of the Veil is affected. This result implies that compo- 
nent B is located closer to the ionized gas than compo- 
nent A, and for the first time reveals the relative locations 
of the two velocity components that constitute the Veil. 
As noted in Sect. 17. ll these two H 1 components are large- 
scale features covering both M42 and M43, and exhibit- 
ing a common velocity gradient, indicating that they are 
related. The fact that they both display a high and sim- 
ilar l ine-of-sight magnetic field strength (iTroland et al. 
1989) reinforces this conclusion. As shown by I Abel et al. 
(2004), the magnetic energy exceeds the energy in tur- 
bulent motions at least in compone nt A, where it con- 
stitutes the dominant pressure term. iFerlandl (|2009f ) has 
suggested that this situation is the result of compression 
by the expanding H 11 region and radiation pressure from 
the ONC. 

The Hi column densities in components A and B, 
based on Lya absorption measurements, a re approxi- 
mate ly 1.6 and 3.2 • 10 21 cur 2 respectively (|Abel et al.l 
2006). With the best density estimates of approximately 
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Fig. 22. — Schematic geometry of M42 and the neutral Veil, illustrating the line-of-sight depths (approximately but not exactly to scale) of 
the various components. Orange indicates ionized gas, green indicates molecular gas and blue indicates atomic gas. The outflow originating 
in Orion-S, which powers the kinematics of HH 202 and H I component F is schematically indicated by a set of thin lines. 



10 2 - 5 and 10 3 - 4 cm- 3 for A and B (jAbel et al.l 12001) . 
their thicknesses are approximately 1.3 (A) and 0.5 pc 
(B). The resulting geometry is schematically shown in 
Fig. The lateral extent of the Veil is at least 1.5 pc 
(its total extent covering M42 and M43), but probably 
significantly larger, since it represents the atomic en- 
velope of the background molecular cloud (vdWG89), 
which extends over at least 10 pc, approximately the 
lengt h of integra l -shap e d filament of molecular ga s and 
dust lBallv et all (fl987t) : [Johnstone fc Ballvl (fl999h . The 
total depth of the ionized cavity is about 0.7 pc, which 
is small compared both to the lateral extent of the Veil 
and to its depth of approximately 1.8 pc. These results 
therefore demonstrate the blister-type nature of the M42 
H II region, emphasized by main IF behind the Trapez- 
ium stars, and also reveal its relation to the envelope of 
the background molecular cloud, represented by the Veil. 

9. CONCLUSIONS 

High resolution H I observations of the Orion Nebula 
reveal in detail the feedback of the H II region and the 
ONC on their neutral environment. This feedback takes 
place through both radiative and mechanical interaction. 
Furthermore, our data provide new information on the 
geometry of the complex. 

Radiative feedback manifests itself principally by the 
non-ionizing UV radiation from the ONC, which creates 
a prominent PDR outside the IF at the interface between 
M42 and OMC-1. Prominent Hi emission is detected 
from the Orion Bar PDR, southeast of the Bright Bar 
IF. This result enables for the first time a study of the 
abundance profile of atomic hydrogen and its relation 
to other tracers in an edge-on PDR. We find that the 
H I emission arises from a region with a temperature of 
~ 540 K, with ~ 5 — 10% of the hydrogen nuclei in the 
interclump region in the PDR in the form of Hi. Most 
of the hydrogen in the region probed therefore remains 



molecular, in agreement with theoretical PDR models. 
This result shows that dense PDRs such as the Orion 
Bar PDR have a low H I production, even if the imping- 
ing UV radiation field is very strong. The Hi column 
density peak matches the peak of UV-excited vibrational 
H2 line emission. This result confirms the two-step H2 
photodissociation process through UV line absorption in 
the Lyman and Werner bands. 

Mechanical feedback results from the interaction of 
ionized and molecular outflow systems with the neutral 
environment of the H II region. The H 1 data provides 
direct evidence for interaction of the HH object HH 202 
and the flow system respons i ble fo r the optical South- 
west Shock of lHennev etaTl (|2007D with the Veil. The 
interaction of these flow systems with the Veil causes ve- 
locity shifts in the absorbing Hi. The evidence for direct 
interaction between HH 202 and the neutral Veil implies 
that the Veil must be located quite close to the Orion-S 
source of outflows. Using the known space trajectory of 
HH 202, we find that the surface of Veil facing the Orion 
Nebula is located at most 0.4 pc in front of Orion-S. Of 
the two velocity components A and B that constitute the 
Veil, we find that component B lies closer to the ionized 
nebula. 

The total depth of the ionized cavity is approximately 
0.7 pc. This depth is small compared to both the lateral 
extent and the depth of Veil. These dimensions confirm 
the nature of the M42 H11 region as a thin blister, and 
simultaneously reveal the relation of the ionized region to 
the background molecular cloud and its neutral envelope. 
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APPENDIX 

RELATION BETWEEN Hi COLUMN DENSITY AND 21 CM LINE OPTICAL DEPTH AND 

BRIGHTNESS TEMPERATURE 

The Hi column density is related to the velocity-integrated optical depth of the 21 cm line by 

^■'-scsifc/*'*' (A1) 

where A ul = 21.10612 cm is the wavelength of the transition, A ul = 2.88426 • 10~ 15 s" 1 (|Gouldlll994l) is its Einstein A 
value, T s is the spin temperature (i.e., the excitation temperature of the two spin levels), h is Planck's constant, k is 
the Boltzmann constant, and c is the speed of light. This expression makes use of the approximation exp(— E n \/kT s ) rs 
1 — E u \/kT s , where E u \ is the energy difference of the two spin levels. This approximation is valid to high accuracy, 
since E u \/k = 0.06816 K and T s 3> E u \/k under all practical conditions. Inserting numbers, Eq. (|A1|) can be written 
as 

= 1.81267. 10 18 § [r(v)d-^. (A2) 
cm z K J km s 1 

For optically thin (t(v) <C 1) Hi emission, exp [— t(v)] can be expanded to first order and the equation of transfer 
Eq. (TTJ) reduces to 

T h (v) = t(v)[T s -T c - T back (v)]. (A3) 
If T s ^> T c + Tback(v) in Eq. (|A3|) . the Hi column density as given by Eq. (|A1|) can be expressed as 

w < Hi >=dSSk/ T *>< fo - < A4 > 

which results in 

^ = 1.81267-10-4^. (A5) 
cm J K km s 1 



REFERENCES 



Abel, N. P., Brogan, C. L., Ferland, G. J., O'Dell, C. R., Shaw, 

G., & Troland, T. H. 2004, ApJ, 609, 247 
Abel, N. P., Ferland, G. J., O'Dell, C. R., Shaw, G., & Troland, 

T. H. 2006, ApJ, 644, 344 
Allen, D. A., & Burton, M. G. 1993, Nature, 363, 54 
Allers, K. N., Jaffe, D. T., Lacy, J. H., Draine, B. T., & Richter, 

M. J. 2005, ApJ, 630, 368 
Ascasibar, Y., Obreja, A. C, & Diaz, A. I. 2011, MNRAS. 1070 
Balick, B., Gammon, R. H., & Hjellming, R. M. 1974, PASP, 86, 

616 

Bally, J. 2008, in Handbook of Star Forming Regions, Volume I, 

ed. Reipurth, B. (ASP), 459-+ 
Bally, J., Cunningham, N. J., Moeckel, N., Burton, M. G., Smith, 

N., Frank, A.. & Nordlund, A. 2011. ApJ, 727, 113 
Bally, J., Langer, W. D., Stark, A. A., & Wilson, R. W. 1987, 

ApJ, 312, L45 

Bally, J., O'Dell, C. R., & McCaughrean, M. J. 2000, AJ, 119, 
2919 

Blaauw, A. 1964, ARA&A, 2, 213 

Black, J. H., & Van Dishoeck, E. F. 1987, ApJ, 322, 412 
Blouin, D., McCutcheon, W. H., Dewdney, P. E., Roger, R. S., 

Purton, C. R., Kester, D., & Bontekoe, T. R. 1997, MNRAS, 

287, 455 
Boisse, P. 1990, A&A, 228, 483 

Brcgman, J., Larson, K., Rank, D., & Temi, P. 1994, ApJ, 423, 
326 

Buckle, J. V., et al. 2012, MNRAS, 422, 521 
Burton, M. G., Hollenbach, D. J., & Tielens, A. G. G. M. 1990, 
ApJ, 365, 620 

Castets, A., Duvert, G., Dutrey, A., Bally, J., Langer, W. D., & 

Wilson, R. W. 1990, A&A, 234, 469 
Chen, H., Bally, J., O'dell, C. R., McCaughrean, M. J., 

Thompson, R. L., Rieke, M., Schneider, G., & Young, E. T. 

1998, ApJ, 492, L173 
Chernin, L. M., & Wright, M. C. H. 1996, ApJ, 467, 676 
Chromcy, F. R., Elmegreen, B. G., & Elmegreen, D. M. 1989, AJ, 

98, 2203 

Clark, B. G. 1965, ApJ, 142, 1398 



Clark, B. G., Radhakrishnan, V., & Wilson, R. W. 1962, ApJ, 
135, 151 

Deharveng, L. 1973, A&A, 29, 341 

Dewdney, P. E., & Roger, R. S. 1982, ApJ, 255, 564 

— . 1986, ApJ, 307, 275 

Doi, T., O'Dell, C. R., & Hartigan, P. 2002, AJ, 124, 445 

— . 2004, AJ, 127, 3456 

Draine, B. T. 1978, ApJS, 36, 595 

Draine, B. T., & Bertoldi, F. 1996, ApJ, 468, 269 

Ferland, G. J. 2009, in IAU Symposium, Vol. 259, IAU 

Symposium, 25—34 
Field, G. B., Somerville, W. B., & Dressier, K. 1966, ARA&A, 4, 

207 

Garcia-Dfaz, M. T., & Henney, W. J. 2007, AJ, 133, 952 
Garcia-Dfaz, M. T., Henney, W. J., Lopez, J. A., & Doi, T. 2008, 

Rev. Mex. Astron. Astr., 44, 181 
Gautier, III, T. N., Fink, U., Larson, H. P., & Treffers, R. R. 

1976, ApJ. 207, L129 
Genzel, R., Reid, M. J., Moran, J. M., & Downes, D. 1981, ApJ, 

244, 884 

Genzel, R., & Stutzki, J. 1989, ARA&A, 27, 41 
Gezari, D. Y., Backman, D. E., & Werner, M. W. 1998, ApJ, 509, 
283 

Goddi, C, Humphreys, E. M. L., Greenhill, L. J., Chandler, 

C. J., & Matthews, L. D. 2011, ApJ, 728, 15 
Goss, W. M., & Shaver, P. A. 1970, Australian Journal of Physics 

Astrophysical Supplement, 14, 1 
Gould, R. J. 1994, ApJ, 423, 522 

Graham, M. F., Meaburn, J., & Redman, M. P. 2003, MNRAS, 
343 419 

Green, D. A. 1991, MNRAS, 253, 350 

Giidel, M., Briggs, K. R., Montmerle, T., Audard, M., Rebull, L., 

& Skinner, S. L. 2008, Science, 319, 309 
Hayashi, M., Hasegawa, T., Gatley, I., Garden, R., & Kaifu, N. 

1985, MNRAS, 215, 31P 
Henney, W. J., & O'Dell, C. R. 1999, AJ, 118, 2350 
Henney, W. J., O'Dell, C. R., Meaburn, J., Garrington, S. T., & 

Lopez, J. A. 2002, ApJ, 566, 315 



28 



Van der Werf et al. 



Henney, W. J., O'Dcll, C. R., Zapata, L. A., Garci'a-Di'az, M. T., 

Rodriguez, L. F., & Robberto, M. 2007, AJ, 133, 2192 
Herrmann, F., Madden, S. C, Nikola, T., Poglitsch, A., 

Timmcrmann, R., Geis, N., Townes, C. H., & Stacey, G. J. 

1997, ApJ, 481, 343 
Heyer, M. H., Morgan, J., Schloerb, F. P., Snell, R. L., & 

Goldsmith, P. F. 1992, ApJ, 395, L99 
Hillenbrand, L. A., & Hartmann, L. W. 1998, ApJ, 492, 540 
Hirota, T., et al. 2007, PASJ, 59, 897 

Hogcrhcijde, M. R., Jansen, D. J., & van Dishoeck, E. F. 1995, 
A&A, 294, 792 

Hollenbach, D. J., & Ticlcns, A. G. G. M. 1999, Reviews of 

Modern Physics, 71, 173 
Huygens, C. 1659, Systcma Satvrnivm (The Hague: A. Vlacq) 
Jaffe, D. T., & Pankonin, V. 1978, ApJ, 226, 869 
Jansen, D. J., Spaans, M., Hogcrhcijde, M. R., & van Dishoeck, 

E. F. 1995, A&A, 303, 541 
Johnston, K. J., Palmer, P., Wilson, T. L., & Bicging, J. H. 1983, 

ApJ, 271, L89 
Johnstone, D., & Bally, J. 1999, ApJ, 510, L49 
Johnstone, D., Boonman, A. M. S., & van Dishoeck, E. F. 2003, 

A&A, 412, 157 

Joncas, G., Durand, D., & Roger, R. S. 1992, ApJ, 387, 591 
Kaler, J. B. 1967, ApJ, 148, 925 

Kassis, M., Adams, J. D., Campbell, M. F., Deutsch, L. K., Hora, 
J. L., Jackson, J. M., & Tollestrup, E. V. 2006, ApJ, 637, 823 

Kaufman, M. J., Wolfire, M. G., Hollenbach, D. J., & Luhman, 
M. L. 1999, ApJ, 527, 795 

Klcinmann, D. E., & Low, F. J. 1967, ApJ, 149, L1+ 

Kraus, S., et al. 2007, A&A, 466, 649 

Kutner, M. L., Tucker, K. D., Chin, G., & Thaddeus, P. 1977, 

ApJ, 215, 521 
Kwan, J., & Scovillc, N. 1976, ApJ, 210, L39 

Le Petit, F., Nehme, C, Le Bourlot, J., & Roueff, E. 2006, ApJS, 
164, 506 

Lis, D. C, & Schilke, P. 2003, ApJ, 597, L145 
Lockhart, I. A., & Goss, W. M. 1978, A&A, 67, 355 
Loren, R. B. 1979, ApJ, 234, L207 

Luhman, K. L., Engelbracht, C. W., & Luhman, M. L. 1998, ApJ, 
499, 799 

Lykins, M. L., Abel, N., & Troland, T. H. 2010, ApJ, 713, 950 
Maddalcna, R. J., Morris, M., Moscowitz, J., & Thaddeus, P. 

1986, ApJ, 303, 375 
Mangum, J. G., Wootten, A., & Plambcck, R. L. 1993, ApJ, 409, 

282 

Marconi, A., Testi, L., Natta, A., & Walmsley, C. M. 1998, A&A, 
330, 696 

Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 
Meijerink, R., Spaans, M., & Israel, F. P. 2007, A&A, 461, 793 
Mcixner, M., & Tielens, A. G. G. M. 1993, ApJ, 405, 216 
Mcntcn, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, 
A&A, 474, 515 

Mcsa-Dclgado, A., Esteban, C, Garcla-Rojas, J., Luridiana, V., 
Bautista, M., Rodriguez, M., Lopez-Martm, L., & Peimbcrt, M. 
2009, MNRAS, 395, 855 
Mezger, P. G., Zylka, R., & Wink, J. E. 1990, A&A, 228, 95 
Mills, B. Y., & Shaver, P. A. 1968, Australian Journal of Physics, 
21, 95 

Muench, A., Getman, K., Hillenbrand, L., & Preibisch, T. 2008, 

in Handbook of Star Forming Regions, Volume I, ed. Rcipurth, 

B. (ASP), 483-+ 
Muller, C. A. 1959, in IAU Symposium, Vol. 9, URSI Symp. 1: 

Paris Symposium on Radio Astronomy, ed. R. N. Braccwell, 360 
Nadcau, D., & Geballe, T. R. 1979, ApJ, 230, L169 
Natta, A., Walmsley, C. M., & Tielens, A. G. G. M. 1994, ApJ, 

428, 209 

O'Dell, C. R. 2001, ARA&A, 39, 99 

O'Dell, C. R., & Harris, J. A. 2010, AJ, 140, 985 

O'Dcll, C. R., Hartigan, P., Lane, W. M., Wong, S. K., Burton, 

M. G., Raymond, J., & Axon, D. J. 1997, AJ, 114, 730 
O'Dell, C. R., & Henney, W. J. 2008, AJ, 136, 1566 
O'Dell, C. R., Henney, W. J., Abel, N. P., Ferland, G. J., & 

Arthur, S. J. 2009, AJ, 137, 367 
O'Dell, C. R., Muench, A., Smith, N., & Zapata, L. 2008, in 

Handbook of Star Forming Regions, Volume I, ed. Reipurth, B. 

(ASP), 544-+ 

O'Dcll, C. R., Valk, J. H., Wen, Z., & Meyer, D. M. 1993, ApJ, 
403, 678 

O'Dcll, C. R., Walter, D. K., & Dufour, R. J. 1992, ApJ, 399, L67 
O'Dcll, C. R., & Wen, Z. 1992, ApJ, 387, 229 
O'Dell, C. R., & Wong, K. 1996, AJ, 111, 846 
O'Dcll, C. R., & Yusef-Zadch, F. 2000, AJ, 120, 382 
Oort, J. H., & Spitzer, Jr., L. 1955, ApJ, 121, 6 



Pellegrini, E. W., Baldwin, J. A., Ferland, G. J., Shaw, G., & 

Heathcote, S. 2009, ApJ, 693, 285 
Plume, R., et al. 2000, ApJ, 539, L133 

Radhakrishnan, V., Brooks, J. W., Goss, W. M., Murray, J. D., 

& Schwarz, U. J. 1972, ApJS, 24, 1 
Rodriguez-Franco, A., Martin-Pintado, J., & Wilson, T. L. 1999a, 

A&A, 351, 1103 
— . 1999b, A&A, 344, L57 

Roger, R. S., & Irwin, J. A. 1982, ApJ, 256, 127 

Roger, R. S., McCutcheon, W. H., Purton, C. R., & Dewdncy, 

P. E. 2004, A&A, 425, 553 
Roger, R. S., & Pedlar, A. 1981, A&A, 94, 238 
Rubin, R. H., Simpson, J. P., O'Dell, C. R., McNabb, I. A., 

Colgan, S. W. J., Zhuge, S. Y., Ferland, G. J., & Hidalgo, S. A. 

2011, MNRAS, 410, 1320 
Sakamoto, S., Hayashi, M., Hasegawa, T., Handa, T., & Oka, T. 

1994, ApJ, 425, 641 
Salas, L., et al. 1999, ApJ, 511, 822 

Schultz, A. S. B., Colgan, S. W. J., Erickson, E. F., Kaufman, 

M. J., Hollenbach, D. J., O'dell, C. R., Young, E. T., & Chen, 

H. 1999, ApJ, 511, 282 
Shaw, G., Ferland, G. J., Henney, W. J., Stancil, P. C, Abel, 

N. P., Pellegrini, E. W., Baldwin, J. A., & van Hoof, P. A. M. 

2009, ApJ, 701, 677 
Shuping, R. Y., & Snow, T. P. 1997, ApJ, 480, 272 
Simon, R., Stutzki, J., Sternberg, A., & Winnewisser, G. 1997, 

A&A, 327, L9 
Simon-Diaz, S. 2010, A&A, 510, A22+ 

Simon-Diaz, S., Hcrrero, A., Esteban, C, & Najarro, F. 2006, 
A&A, 448, 351 

Stacey, G. J., Jaffe, D. T., Geis, N., Grenzel, R., Harris, A. I., 

Poglitsch, A., Stutzki, J., & Townes, C. H. 1993, ApJ, 404, 219 
Stecher, T. P., & Williams, D. A. 1967, ApJ, 149, L29+ 
Sternberg, A. 1988, ApJ, 332, 400 
Sternberg, A., & Dalgarno, A. 1989, ApJ, 338, 197 
Stolovy, S. R., et al. 1998, ApJ, 492, L151+ 

Subrahmanyan, R., Goss, W. M., & Malin, D. F. 2001, AJ, 121, 
399 

Takami, M., et al. 2002, ApJ, 566, 910 

Taubcr, J. A., Ticlcns, A. G. G. M., Mcixner, M., & Foldsmith, 

P. F. 1994, ApJ, 422, 136 
Tielens, A. G. G. M., & Hollenbach, D. 1985a, ApJ, 291, 747 
— . 1985b, ApJ, 291, 722 

Tielens, A. G. G. M., Meixner, M. M., Van der Werf, P. P., 
Bregman, J., Tauber, J. A., Stutzki, J., & Rank, D. 1993, 
Science, 262, 86 

Troland, T. H., Heiles, C, & Goss, W. M. 1989, ApJ, 337, 342 
Van der Werf, P. P., & Goss, W. M. 1989, A&A, 224, 209 
— . 1990, ApJ, 364, 157 

Van der Werf, P. P., Stutzki, J., Sternberg, A., & Krabbe, A. 

1996, A&A, 313, 633 

Van der Wiel, M. H. D., van der Tak, F. F. S., Ossenkopf, V., 
Spaans, M., Roberts, H., Fuller, G. A., & Plume, R. 2009, 
A&A, 498, 161 

Van Gorkom, J. H., & Ekcrs, R. D. 1989, in Astronomical Society 
of the Pacific Conference Series, Vol. 6, Synthesis Imaging in 
Radio Astronomy, ed. R. A. Perley, F. R. Schwab, & 
A. H. Bridle, 341-+ 

Walmsley, C. M., Natta, A., Oliva, E., & Testi, L. 2000, A&A, 
364, 301 

Weigelt, G., Balega, Y., Preibisch, T., Schertl, D., Scholler, M., & 

Zinnecker, H. 1999, A&A, 347, L15 
Wen, Z., & O'Dell, C. R. 1995, ApJ, 438, 784 
Wilson, B. A., Dame, T. M., Masheder, M. R. W., & Thaddeus, 

P. 2005, A&A, 430, 523 
Wilson, T. L., Filges, L., Codella, C, Reich, W., & Reich, P. 

1997, A&A, 327, 1177 

Wyrowski, F., Schilke, P., Hofncr, P., & Walmsley, C. M. 1997, 

ApJ, 487, L171+ 
Young Owl, R. C, Meixner, M. M., Wolfire, M., Tielens, 

A. G. G. M., & Tauber, J. 2000, ApJ, 540, 886 
Yusef-Zadch, F. 1990, ApJ, 361, L19 

Zapata, L. A., Ho, P. T. P., Rodriguez, L. F., Schilke, P., & 

Kurtz, S. 2007, A&A, 471, L59 
Zapata, L. A., Loinard, L., Schmid-Burgk, J., Rodriguez, L. F., 

Ho, P. T. P., & Patel, N. A. 2011a, ApJ, 726, L12+ 
Zapata, L. A., Schmid-Burgk, J., Ho, P. T. P., Rodriguez, L. F., 

& Menten, K. M. 2009, ApJ, 704, L45 
Zapata, L. A., Schmid-Burgk, J., & Menten, K. M. 2011b, A&A, 

529, A24+ 
Zuckerman, B. 1973, ApJ, 183, 863 

Zuckerman, B., Kuiper, T. B. H., & Rodriguez Kuipcr, E. N. 
1976, ApJ, 209, L137 



